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Abstract 

This  research  analyzes  material  characterization  measurements  where  axially 
symmetric  air  gaps  exist  between  the  sample  material  and  the  inner  or  outer  conductor 
of  a  coaxial  test  fixture.  Higher  order  fields  are  excited  by  the  air  gap  and  are  not 
accounted  for  in  generally  used  algorithms  for  determining  the  material  permittivity 
and  permeability.  The  result  is  error  in  the  material  characterization  measurement. 

This  research  defines  the  fields  within  the  material  and  air  gap,  assuming  them 
to  be  axially  symmetric.  The  fields  are  then  used  via  the  modal  method  to  calcu¬ 
late  theoretical  scattering  parameters  as  a  function  of  permittivity,  permeability,  and 
frequency.  A  complex,  two-dimensional  Newton  root  search  then  iterates  the  per¬ 
mittivity  and  permeability  for  a  given  frequency  minimizing  the  difference  between 
the  calculated  scattering  parameters  and  the  measured  scattering  parameters.  In  this 
manner,  the  root  corresponds  to  the  permittivity  and  permeability  of  the  sample. 

The  modal  method  provides  accurate  results  for  non-magnetic  material  mea¬ 
surements  when  the  material  sample  fills  only  30%  of  the  radial  distance  between  the 
inner  and  outer  conductor  of  the  coaxial  line.  Due  to  the  concentration  of  the  electric 
field  at  the  inner  conductor,  accurate  results  were  achieved  with  a  29.25  mil  material 
layer  (400  mil  air  gap)  on  the  inner  conductor.  It  is  shown  that  the  modal  method 
result  converges  to  the  material  properties  by  using  10  modes. 

The  modal  method  provides  good  results  for  high-dielectric  constant  magnetic 
material.  Results  for  an  outer  gap  scenario  were  more  accurate  than  inner  gap  results. 
The  modal  method  appears  to  converge  with  very  little  error  at  approximately  20 
modes,  but  diverges  for  greater  than  20  modes.  This  behavior  is  not  expected  and 
is  attributed  to  a  failure  of  computer  calculated  Bessel  functions  with  large  real  and 
small  complex  arguments. 
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The  sensitivity  of  results  to  uncertainty  in  sample  radial  thickness  and  length 
measurements  were  analyzed.  Uncertainty  in  the  radial  thickness  was  the  dominate 
source  of  error,  particularly  when  the  air  gap  is  large  or  the  material  is  heavily  loaded. 

Two  factors  significantly  impacted  the  performance  of  the  modal  method.  First, 
the  presence  of  an  air  gap  effects  the  physical  support  to  the  center  conductor.  For 
the  air  gap  to  remain  axially  symmetric  (a  fundamental  assumption  of  this  research), 
the  center  conductor  must  remain  along  the  axis  of  the  test  fixture  outer  conductor. 
Second,  accurately  determining  the  z-directed  wavenumber  in  the  material  and  air 
gap  region  is  critical  and  requires  much  care,  particularly  when  measuring  heavily 
loaded  materials. 
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AIR  GAP  ERROR  COMPENSATION  FOR 


COAXIAL  TRANSMISSION  LINE  METHOD  OF 
ELECTROMAGNETIC  MATERIAL  CHARACTERIZATION 

I.  Introduction 

A  material’s  response  to  the  application  of  an  electromagnetic  field  is  quantified 
by  the  material’s  complex  permittivity  and  permeability.  The  real  and  imaginary 
parts  can  be  viewed  as  the  susceptibility  to  polarization  and  attenuation  (or  loss) 
respectively.  The  permittivity  is  in  response  to  the  electric  held  and  the  permeability 
to  the  magnetic  held. 

A  material’s  permittivity  and  permeability  are  exploited  for  a  wide  variety  of 
antenna,  electromagnetic  interference,  and  low  signature  applications.  Accurate  char¬ 
acterization  of  a  material’s  electromagnetic  properties  is  therefore  a  necessary  and  cru¬ 
cial  step  in  the  research,  development,  and  design  of  new  material  mixture  formulas 
and  in  the  verification  of  manufacturing  processes. 

1 . 1  Material  Characterization 

The  Transmission/Rehection  (T/R)  method  of  material  characterization  deter¬ 
mines  the  complex  permittivity  and  permeability  of  a  linear,  homogeneous,  isotropic 
material  by  measuring  the  material’s  scattering  parameters  [4],  A  sample  of  the  ma¬ 
terial  is  mounted  internal  to  a  test  hxture  and  an  attached  network  analyzer  measures 
the  scattering  parameters  necessary  to  determine  the  material’s  electromagnetic  prop¬ 
erties.  The  well-known  Nicolson- Ross- Weir  (NRW)  algorithm  extracts  the  material 
parameters  from  the  measured  scattering  parameters  [12,17].  The  NRW  algorithm 
assumes  the  unknown  material  completely  fills  the  cross  sectional  area  of  the  test  hx¬ 
ture.  This  assumption,  in  conjunction  with  the  measurement  frequency  band,  restricts 
the  helds  within  the  test  hxture  to  a  single  propagating  mode. 
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Discontinuities  may  unintentionally  exist  within  the  test  fixture  as  a  result  of 
imprecise  machining  of  the  sample  or  other  defect.  Deliberate  discontinuities  are  also 
possible.  For  example,  high  temperature  measurements  are  prone  to  gaps  due  to  the 
different  thermal  expansion  rates  for  the  material  and  conductor. 

Discontinuities  in  the  test  fixture  will  excite  high-order  evanescent  or  higher- 
order  propagating  modes.  Due  to  the  conservation  of  energy  and  mode  orthogonality, 
the  incident  field  power  will  couple  into  both  classes  of  higher  order  modes.  Detection 
circuitry  within  the  network  analyzer  is  typically  too  far  from  the  forward  face  of 
the  material  under  test  to  detect  power  coupled  into  higher-order  evanescent  waves. 
Additionally,  the  test  setup  (network  analyzer  and  test  fixture)  is  normally  designed 
to  excite  and  detect  only  one  field,  the  dominant  mode.  If  higher-order  propagat¬ 
ing  modes  are  excited,  their  power  will  not  be  effectively  identified  by  the  detector 
element. 

1.2  Accounting  for  Discontinuity  Error 

The  power  loss  to  higher-order  modes  is  a  source  of  measurement  error  in  the 
NRW  algorithm.  This  research  seeks  to  account  for  the  error  induced  when  higher- 
order  modes  are  excited  in  order  to  accurately  determine  the  material  parameters  in 
such  cases. 

The  test  fixture  for  this  research  is  the  coaxial  line  shown  in  figure  1.1.  The 
discontinuities  considered  in  this  research  are  axially  symmetric  air  gaps  between  the 
outer  radius  of  the  sample  or  the  inner  radius  of  the  outer  conductor. 
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II.  Discontinuity  Error  in  Coaxial  Line  Material 

Measurements 

The  T /R  method  of  material  characterization  is  widely  used  due  to  its  simplicity  [4] 
and  is  the  measurement  method  for  this  research.  The  method  requires  a  single 
measurement  device,  typically  a  network  analyzer  (NWA),  and  a  means  to  mount  the 
sample  for  testing,  typically  a  waveguide  test  fixture.  This  research  uses  a  coaxial  line 
as  the  test  fixture.  The  system  is  connected  as  shown  in  figure  2.1. 

2.1  Why  a  Coaxial  Test  Fixture? 

Measurement  frequency  band  is  a  primary  concern  when  selecting  a  test  fixture. 
A  two  conductor,  coaxial  transmission  line’s  ability  to  propagate  a  TEM  wave  makes 
it  an  ideal  fixture  for  tests  that  require  relatively  broadband,  low  frequency  measure¬ 
ments.  Additionally,  the  typically  small  size  of  a  coaxial  line  ensures  a  reasonable 
material  sample  size. 

Each  field  configuration  (mode)  within  a  waveguide  has  a  lower  frequency  bound 
known  as  the  cutoff  frequency.  The  propagating  mode  with  the  lowest  frequency 
bound  is  the  dominant  mode  for  a  given  waveguide.  The  cutoff  frequency  is  de¬ 
termined  from  waveguide  boundary  conditions  and,  in  the  case  of  rectangular  and 
cylindrical  single-conductor  waveguide,  is  inversely  proportional  to  the  dimensions 
of  the  waveguide  therefore  requiring  very  large  dimensions  to  measure  low  frequen¬ 
cies  [3].  On  the  other  hand,  a  multi-conductor  waveguide  is  unique  in  its  ability  to 
propagate  a  TEM  wave.  The  cutoff  frequency  for  a  TEM  wave  is  0  Hz  (see  appendix 
C)  making  it  the  dominant  mode  in  a  multi-conductor  waveguide  such  as  a  coaxial 
line.  This  fact  permits  low  frequency  material  measurements. 

The  dimensions  of  the  coaxial  test  fixture  dictate  the  measurement  frequency 
bandwidth.  The  coaxial  test  fixture  for  this  research  has  an  inner  conductor  of  radius 
a  =  0.331  inches  and  outer  conductor  radius  of  b  =  0.76025  inches.  The  bandwidth 
of  the  measurement  is  typically  kept  at  ~  90%  of  the  dominant  mode  bandwidth.  As 


4 


Figure  2.1:  Test  setup  with  an  exploded  view  of  the  test  fixture  and  mounted 

sample.  The  forward  scattering  parameters,  Sn  and  S21,  are  measured  at  the  NWA 
ports  1  and  2  respectively. 
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seen  from  table  2.1,  the  first  higher  order  mode  in  the  coaxial  test  fixture  is  the  TEn 
with  a  cutoff  frequency  of  3.5  GHz  [10].  This  sets  the  upper  frequency  bound. 


Although  the  cutoff  of  a  TEM  is  a  0  Hz,  the  lower  frequency  bound  of  the  NWA 
is  50  MHz.  Therefore,  the  measurement  frequency  bandwidth  for  this  research  is 
50  Mhz  to  3.05  GHz  or  m  85%  of  the  available  TEM  bandwidth. 

The  frequency  bounds  must  consider  the  evanescent  fields  excited  by  the  transi¬ 
tion  from  the  NWA  cables  to  the  coaxial  test  fixture.  These  fields  may  not  propagate, 
but  they  may  not  decay  rapidly  either.  For  example,  a  50  Mhz  TEM  wave  in  the 
freespace  hlled  region  of  the  coaxial  line  in  this  research  will  lose  63%  of  its  field 
strength  in  approximately  0.94  m.  If  a  sample  is  placed  to  close  to  the  end  of  the 
coaxial  test  fixture,  there  will  be  error  due  to  the  higher  order  evanescent  modes  ex¬ 
cited  at  the  transition.  This  research  assumes  a  single  mode  excitation,  meaning  the 
sample  is  far  enough  away  from  the  end  of  the  coaxial  test  fixture  with  respect  to  the 
given  measurement  frequency  band. 

The  coaxial  line  dimensions  also  define  the  amount  of  sample  needed  and  the  im¬ 
pedance  of  the  test  fixture.  The  coaxial  line  in  this  research  requires  only  a  1.5205  in. 
by  1.5205  in.  sample  to  fill  the  cross  section.  Additionally,  the  coaxial  line  dimensions 
produce  an  impedance 


necessary  to  match  the  characteristic  impedance  of  the  NWA  cables. 

2.2  NRW  Algorithm  and  Discontinuities 

The  network  analyzer  excites  a  field,  typically  of  the  dominant  mode,  that  trav¬ 
els  through  the  waveguide  until  it  strikes  the  mounted  sample.  A  portion  of  the 
incident  field  is  reflected  and  a  portion  transmitted  through  the  material.  The  mag¬ 
nitude  and  phase  of  the  reflected  and  transmitted  field  is  measured  by  the  network 
analyzer  and  reported  as  complex  scattering  parameters  Sn  and  S21  (see  figure  2.1). 
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Tabic  2.1:  The  first  10  modes  in  a  coaxial  line  with  inner  conductor  of  radius 

0.331  inches  and  outer  conductor  radius  of  0.76025  inches.  The  TE02  and  TMn  are 
expected  to  have  the  same  cutoff  frequency  [10]. 


Freespace  Filled  Coax 

Mode 

Cutoff  Frequency  [GHz] 

7 

kp  kc 

TEM 

0.000 

0.0+41. 9j 

O.O+O.Oj 

TEn 

3.514 

60.6 

73.6 

TE21 

6.888 

138.1 

144.4 

TE3i 

10.041 

206.2 

210.4 

TE4i 

12.986 

268.9 

272.2 

TM01 

13.633 

282.6 

285.7 

TE02 

14.088 

292.3 

295.3 

TMn 

14.088 

292.3 

295.3 

te12 

14.652 

304.2 

307.1 

tm21 

15.364 

319.3 

322.0 

Physically,  Sn  and  S21  are  the  ratio  of  the  reflected  field  amplitude  and  transmitted 
field  amplitude,  respectively,  to  the  incident  field  amplitude  under  matched  condi¬ 
tions.  The  NRW  algorithm  provides  analytical  expressions  for  the  permittivity  and 
permeability  of  the  sample  as  a  function  of  the  measured  scattering  parameters1,  Sn 
and  S21  [4, 12, 17]. 

A  calibrated  system  is  critical  for  obtaining  accurate  measurements.  Mismatches 
at  cable  connections,  waveguide  deformities,  and  coupling  between  measurement  de¬ 
vices  at  ports  1  and  2  are  just  a  few  items  that  will  introduce  error  into  the  mea¬ 
surement.  Calibration  will  significantly  reduce  the  impact  these  defects  have  on  the 
measurements.  This  research  assumes  a  properly  calibrated  system  via  a  Thru,  Re¬ 
flect,  Line  (  TRL)  technique  and  that  measured  data  is  representative  of  the  proper¬ 
ties  of  the  sample,  not  system  imperfections.  The  details  of  system  calibrations  and 
specifically  the  TRL  calibration,  are  thoroughly  discussed  in  [4], 

1Although  the  forward  scattering  parameters  Sn  and  S21  are  discussed  in  the  text,  they  are 
interchangeable  with  the  reverse  scattering  parameters  S22  and  S12  respectively.  Measurement  of 
S22  and  S12  require  exciting  the  incident  field  at  port  2. 
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A  fundamental  assumption  of  the  NRW  method  is  that  the  power  in  the  system 
is  carried  by  a  single  mode  (typically  the  dominant  mode)  and  is  therefore  completely 
measured  by  the  NWA  and  quantified  in  the  Sn  and  S21  measurements.  The  or¬ 
thogonality  of  field  modes  (see  appendix  D)  ensures  that  a  single  mode  excitation 
in  a  uniform  guiding  structure  will  carry  the  power  in  the  system.  If  higher  order 
modes  are  excited  within  the  structure,  power  carried  by  the  initial  excitation  will 
couple  into  the  higher  order  modes,  but  will  typically  evanescence  due  to  the  cutoff 
frequency  of  the  higher  mode  being  greater  than  the  measurement  frequency  band. 
This  represents  a  loss  in  power  and,  consequently,  perturbation  in  the  Sn  and  S21 
measurements.  Not  all  discontinuities  will  excite  higher  order  modes  (as  shown  in 
appendix  H),  but  if  higher  order  modes  are  excited,  they  must  be  accounted  for  to 
produce  accurate  material  parameter  results. 

2.3  Radial  Air  Gap  Discontinuities 

Radial  air  gap  discontinuities  as  shown  in  figure  2.2b  and  2.2c  are  the  primary 
focus  of  this  research.  Samples  are  typically  machined  via  a  lathe  by  experienced 
machinists  and  therefore,  if  discontinuities  exist,  they  are  assumed  to  be  axially  sym¬ 
metric.  Samples  are  machined  to  fill  the  entire  cross-sectional  area  of  the  test  fixture 
(figure  2.2a)  for  the  length  of  the  sample  and  are  press  fit  into  the  coaxial  test  fix¬ 
ture  of  figure  1.1,  minimizing  the  potential  for  discontinuities  and  higher-order  mode 
excitation. 

Repeated  mounting  of  a  single  sample,  imprecise  machining,  or  malleable  mate¬ 
rial  may  unintentionally  introduce  discontinuities  between  the  sample  and  test  fixture 
conductors.  Each  time  a  material  is  mounted,  some  of  the  outer  and  inner  surface  is 
removed.  Optimally,  the  amount  of  material  removed  is  kept  to  a  minimum  by  care¬ 
ful  alignment  of  the  sample  during  mounting.  Precision  machining  is  necessary,  but 
difficult  and  potentially  expensive.  Unintentional  discontinuities  are  usually  small  in 
scale  and  can  be  compensated  for  with  conducting  pastes  between  the  material  and 
conductor  [15],  but  risk  contaminating  the  sample. 


(c)  (d) 

Figure  2.2:  A  sample  that  uniformly  fills  the  cross-sectional  area  of  the  coaxial  test 
fixture,  shown  in  (a),  is  the  ideal  test  scenario.  Machining  imperfections  or  temper¬ 
ature  changes  during  measurement  can  introduce  air  gaps  between  the  outer  coaxial 
conductor  and  the  outer  radius  of  the  sample  (b)  or  between  the  inner  conductor  and 
the  inner  radius  of  the  sample  (c).  A  two  air  gap  generalization  (d)  could  be  used  for 
analysis,  but  physically  would  not  remain  axially  symmetric.  Scenarios  (b)  and  (c) 
are  addressed  in  this  research. 
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HEAT 


Ceramic 
Material 

At  Room  Temperature  At  Measurement  Temperature 


Figure  2.3:  High  temperature  measurements  of  ceramic  materials  may  crack  the 

material  if  the  inner  conductor  expands  faster  than  the  material  as  temperature  in¬ 
creases  (upper).  Machining  in  an  inner  air  gap  (lower  left)  ensures  the  sample  is 
undamaged  at  measurement  temperature,  but  also  introduces  an  outer  air  gap  (lower 
right)  and  consequently  error  in  the  measurement. 

Test  profiles  may  require  that  relatively  large  gaps  be  machined  into  the  material 
prior  to  measurement.  For  example,  measuring  a  ceramic  material  at  high  temper¬ 
ature  requires  the  sample  be  machined  with  an  inner  air  gap  at  room  temperature 
(figure  2.3).  The  mismatch  in  thermal  expansion  rates  between  the  sample  and  the 
conductor  may  crack  or  otherwise  damage  the  material  at  measurement  temperature. 
If  the  test  profile  calls  for  measurements  at  high  temperatures,  a  sample  may  require 
an  inner  air  gap  at  room  temperature  that  is  filled  in  as  the  coax  center  conductor 
expands  with  increase  in  temperature  (assuming  a  mismatch  in  the  thermal  expan¬ 
sion  coefficients  between  the  coax  metal  and  the  sample).  Obviously,  as  temperature 
increases  the  outer  conductor  is  also  expanding  away  from  the  outer  material  face 
leaving  an  unavoidable  outer  air  gap. 
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The  results  from  measurements  of  two  different  materials  are  presented  as  an 
example  of  the  impact  of  air  gaps.  Both  materials  are  assumed  to  be  a  simple  media 
(i.e.  linear,  isotropic,  and  homogeneous).  The  first  material  is  Rexolite®,  a  non¬ 
magnetic,  lossless  dielectric  material  with  a  relative  permittivity  of  approximately 
2.525  over  the  measurement  frequency  band.  Three  samples  were  used: 

1.  Fully-Filled:  the  sample  fully- fills  the  cross-sectional  area  of  the  coaxial  waveguide 
between  the  inner  and  outer  conductor.  The  sample  was  press-fitted  dry  into 
the  sample  holder  shown  in  figure  1.1.  Measurements  from  this  sample  are  as¬ 
sumed  to  be  representative  of  the  true  properties  of  the  material  given  the  test 
setup  used. 

2.  100  mil  Outer  Gap:  the  sample  was  machined  to  have  a  100  mil  radial  air  gap 
between  the  outer  conductor  and  the  material.  The  sample  was  press-fitted  dry 
onto  the  inner  conductor.  The  air  gap  is  ~  30%  of  the  cross-sectional  area  in 
the  coax  test  fixture  between  the  inner  and  outer  conductor. 

3.  100  mil  Inner  Gap:  the  sample  was  machined  to  have  a  100  mil  radial  air  gap 
between  the  inner  conductor  and  the  material.  The  sample  was  press-fitted  dry 
into  the  outer  conductor.  The  air  gap  is  ~  17%  of  the  cross-sectional  area  in 
the  coax  test  fixture  between  the  inner  and  outer  conductor. 

Each  sample  was  1  in.  in  length. 

The  second  material  is  Magnetic  Radar  Absorbing  Material  (MagRAM),  a  high 
dielectric  constant,  lossy,  magnetic  material.  It  is  designed  for  applications  requiring 
significant  attenuation  of  electric  and  magnetic  fields.  The  permittivity  and  perme¬ 
ability  of  the  MagRAM  vary  with  frequency  and  are  shown  in  subsequent  figures  as 
the  fully  Elled  result.  Three  MagRAM  samples  were  used  and  are  identical  in  de¬ 
scription  to  the  Rexolite®  samples  except  for  their  radial  and  length  measurements. 
One  MagRAM  sample  fully  filled  the  cross-section  of  the  coaxial  test  fixture,  while 
the  other  two  samples  had  9  mil  inner  and  outer  air  gaps.  The  9  mil  outer  air  gap  is 
3%  of  the  cross  sectional  area  of  the  coaxial  test  fixture  between  the  inner  and  outer 
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conductor.  The  9  mil  inner  air  gap  is  1.3%.  The  MagR.AM  sample  lengths  for  the 
fully  filled,  outer  and  inner  gap  are  173.4  mil,  150  mil  and  179  mil  respectively.  The 
MagRAM  samples  were  cut  from  sheets  of  the  material  with  care  taken  to  minimize 
variations  in  the  material  composition  (i.e.  cutting  samples  in  the  same  general  area 
of  the  sheet). 

The  impact  of  air  gaps  on  the  permittivity  measurements  follows  physical  in¬ 
tuition.  The  permittivity  measurements  for  Rexolite®  are  shown  in  figure  2.4.  As 
evident  in  (a),  introducing  an  air  gap  in  the  measurement  is  equivalent  to  “mixing” 
air  with  the  material  to  be  measured.  Therefore,  the  real  permittivity  is  expected 
to  decrease,  approaching  that  of  freespace  (eo,/io)  for  large  gaps.  The  measurements 
validate  this  expectation.  The  imaginary  permittivity  is  expected  to  remain  constant 
(the  material  and  freespace  are  both  very  low  loss).  The  NRW  algorithm  results  do 
not  validate  this  expectation.  The  coupling  of  power  into  higher  order  evanescent 
modes  represents  a  “loss”  in  system  power  as  measured  at  the  NWA.  The  NRW  al¬ 
gorithm  expects  that  the  NWA  measurements  (Sn  and  S21)  are  representative  of  the 
sample  permittivity  and  permeability  and  therefore  misreports  the  sample  as  lossy. 
In  neglecting  the  higher  order  mode  excitation,  the  NRW  algorithm  is  accurately  re¬ 
porting  a  property  of  the  system,  but  not  of  the  material  being  measured.  Similar 
results  are  noted  for  the  MagRAM  material  shown  in  figure  2.5. 

The  location  of  the  air  gap  impacts  the  degree  of  error  in  the  measurement. 
The  electric  field  is  concentrated  around  the  center  conductor  of  the  coaxial  line.  To 
measure  the  effect  of  the  electric  and  magnetic  field  on  a  material  (permittivity  and 
permeability),  the  material  is  optimally  exposed  to  the  field  strengths  of  each  field.  For 
inner  air  gaps,  the  freespace  region  is  exposed  to  the  peak  field  strengths.  Therefore, 
a  greater  error  is  expected  when  an  air  gap  exists  around  the  inner  conductor  versus 
an  equivalent  width  air  gap  around  the  outer  conductor.  The  measurements  shown 
in  figure  2.4  and  2.5  match  this  expectation. 
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Figure  2.4:  Measured  permittivity  is  shown  as  er  =  e'r  —  je”.  Rexolite®  is  a  non¬ 
magnetic,  lossless  dielectric  (er  ~  2.525).  The  real  part  of  the  permittivity  (a)  de¬ 
creases  toward  that  of  freespace  when  an  air  gap  is  introduced.  Additionally,  the 
imaginary  part  of  the  permittivity  (b)  indicates  more  loss  when  higher  order  modes 
are  excited,  but  not  accounted  for.  The  concentration  of  the  electric  held  around  the 
center  gives  rise  to  more  error  if  an  air  gap  exists  around  the  inner  conductor. 


2.4  Air  Gap  Correction 


Frequency  independent  permittivity  and  permeability  corrections  for  coaxial 
line  measurements  where  an  air  gap  exists  in  the  sample  region  are  available  from  the 
National  Institute  of  Standards  and  Technology  (NIST)  [2],  The  corrections  assume 
that  the  fields  in  the  sample  region  are  TEM.  The  sample  with  air  gap  is  then  modelled 
independently  as  a  layered  capacitor  and  a  series  inductor  where  the  singely  directed 
components  of  the  electric,  Ep,  and  magnetic,  fields  are  applied  across  the  lumped 
circuit  components.  As  an  example,  using  the  generalized  case  of  2.2(d),  the  total 
capacitance  of  the  sample  with  air  gaps  is  equivalent  to  a  layered  or  series  capacitor 
and  is  given  by 


a 


total 


a 


+ 


outer  air  gap 


a 


sample 


a 


inner  air  gap 


(2.1) 
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- Fully  Filled  - 9  mil  Outer  Gap  - 9  mil  Inner  Gap 


Frequency  [GHz]  Frequency  [GHz] 

(a)  (b) 


Figure  2.5:  Measured  permittivity  and  permeability  of  MagRAM  are  shown  as 

er  =  e'r  —  je"  and  fir  =  n'r  —  j/i ",  respectively.  Even  small  air  gaps  have  a  significant 
impact  on  the  measurement  of  highly  load  materials.  The  exact  effects  are  as  noted 
before,  a  decrease  in  measured  real  permittivity  and  permeability  toward  that  of 
freespace  and  an  increase  in  imaginary  part  of  both  the  permittivity  and  permeability 
due  to  unaccounted  for  power  coupled  into  higher  order  evanescent  modes. 
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Knowing  that  capacitance  is  the  ratio  of  charge  to  potential,  the  capacitance  of  each 
individual  layer  is  equivalent  to  [6] 


C 


Q 

V 

e  <f  p  ■  pEp  dl 


f  pEp  ■  dl 

pi 


where  surface  S  is  the  cross  sectional  area  of  the  layer,  pi  and  P2  are  the  inner  and 
outer  radius  of  the  layer  respectively,  e  is  the  permittivity  of  the  layer,  L  is  the  length 
of  the  sample,  and  pEp  is  the  p  directed  electric  field  for  a  TEM  excitation  or 
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Solving  C  in  a  radial  layer  produces 
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which  can  be  calculated  for  each  layer  and  substituted  into  (2.1)  to  produce  simple, 
analytical  correction  factors  for  the  real  and  imaginary  parts  of  the  permittivity  cal¬ 
culated  via  the  NRW  algorithm.  An  identical  process  is  accomplished  for  the  series 
inductance  and  results  in  a  correction  factor  for  the  permeability.  The  corrections 
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are  [2] 
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where  e^R  —  j^'mR  and  n'mR  —  jl^mR  are  tlie  relative  permittivity  and  permeability 
calculated  from  measurements,  e'cR  —  je"R  and  n'cR  —  jn"R  are  the  corrected  relative 
permittivity  and  permeability,  and  using  dimensions  from  figures  2.2b  and  2.2c 


The  correction  factors  are  very  simple  to  apply  to  measured  data  and  require  only 
knowledge  of  the  radii  of  the  coaxial  line  and  sample. 

In  1994,  NIST  enlisted  eleven  organizations  to  independently  measure  five  unique, 
low- loss  material  samples  of  known  permittivity  to  compare  each  organizations  method 
of  air  gap  error  correction2.  The  real  part  of  the  permittivity,  er  ',  for  the  five  mate¬ 
rials  was  on  the  order  of  10  with  one  notable  exception  of  er  '  —  50  (Barium  titanate 
ceramic).  Based  on  participant  measurements,  the  samples  had  inner  and  outer  air 
gaps  on  the  order  of  1  mil  so  that  the  air  gap  accounted  for  2.5%  of  the  cross  sectional 
area  of  the  coaxial  line  while  the  sample  filled  the  remaining  97.5%.  Participants  used 
one  of  three  means  to  calculate  the  material  permitttivity:  the  NRW  algorithm  with 
no  corrections,  the  NRW  algorithm  with  the  NIST  correction  and  other  undisclosed 
(assumed  proprietary)  techniques.  The  results,  published  in  [15],  show  most  mea¬ 
surements  agreeing  within  10%,  although  a  great  deal  of  frequency  dependent  error 

2Measurements  were  taken  at  room  temperature  and  the  use  of  conductive  pastes,  mentioned 
previously,  was  forbidden. 
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is  noted  in  the  high  permittivity  material.  The  authors  concluded  that  “accurate 
knowledge  of  the  air-gap  dimensions  is  fundamental  to  the  proper  determination  of 
material  characteristics  using  the  T /R  method.” 

The  NIST  correction  factors  were  applied  to  the  Rexolite®  and  MagRAM  mea¬ 
surements  to  judge  their  effectiveness  on  the  materials  used  for  this  research.  The 
resulting  corrected  measurements  are  shown  in  figures  2.6  and  2.7  respectively. 

Corrected  e'r  values  for  both  the  Rexolite®  and  MagRAM  are  reasonably  ac¬ 
curate.  Figure  2.8  provides  more  detail  of  the  MagRAM  e'r  corrections  at  2  GHz 
and  showing  a  general  trend  toward  the  vicinity  of  the  actual  value  despite  the  large 
discrepancy  in  low  frequency  corrected  values.  In  the  case  of  MagRAM,  fj!r  is  under¬ 
corrected,  but  the  correction  is  toward  the  actual  value  of  fj!r  and  maintains  the  same 
general  slope  over  a  change  in  frequency. 

The  corrected  imaginary  material  parameters  do  not  give  a  good  understanding 
of  the  behavior  of  each  material.  For  Rexolite®,  the  corrected  e"  is  diverging  from  the 
actual  material  parameters.  For  MagRAM,  both  the  corrected  e”  and  /z"  are  diverging 
from  the  actual  value.  These  results  are  expected  since  higher  order  evanescent  fields 
are  not  being  accounted  for  in  either  the  NRW  algorithm  or  the  NIST  correction. 

Herein  lies  the  difficulty  in  applying  the  correction  factor.  For  both  the  low  and 
high  permittivity  materials,  the  amount  of  correction  is  fixed.  Although  a  sensitivity 
analysis  may  lend  confidence  to  the  final  result,  no  physical  insight  is  readily  available 
from  the  correction  factor.  One  of  the  acknowledged  shortcomings  of  the  layered 
capacitor  model  is  the  assumption  of  a  single-propagating  mode  and  that  the  fields 
internal  to  the  material  region  are  TEM.  It  is  already  proven  that  higher  order  modes 
will  be  excited  when  an  air  gap  discontinuity  is  present.  It  is  now  necessary  to  define 
the  fields  within  the  material  region  to  determine  how  close  they  are  to  TEM. 
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- Fully  Filled 


100  mil  Outer  Gap  -  No  Fix 

100  mil  Outer  Gap  -  NIST  Correction 


100  mil  Inner  Gap  -  No  Fix 

100  mil  Inner  Gap  -  NIST  Correction 


Figure  2.6:  Measured  permittivity  is  shown  as  er  =  e'r  —  je".  Rexolite®  is  a  non¬ 
magnetic,  lossless  dielectric  (er  ~  2.525).  The  NIST  correction  performs  very  well  for 
the  e'r  (a),  but  e"  is  over-corrected  (b). 


2.5  Fields  in  Region  with  Sample  and  Air  Gap 

Air  gaps  between  the  sample  and  the  outer  or  inner  conductor  of  the  coaxial  test 
fixture  eliminate  the  possibility  of  a  TEM  held  within  the  material.  Considering  a  z- 
directed  held  as  in  figure  2.9,  the  phase  front  must  propagate  uniformly  through  the  air 
gap  and  material  so  that  helds  at  the  material  to  air  gap  interface  are  continuous,  or, 
mathematically,  independent  of  z.  For  the  l-directed  phase  front  to  remain  continuous 
across  the  radial  boundary 

e-jk°zz  =  e-jk%z  ( 2 ,2) 

must  be  true,  where  kz  and  ksz  are  the  ^-directed  wave  numbers  of  free  space  and  the 
sample  respectively.  The  cutoff  wave  number  of  a  TEM  held  is  kc  =  krh0  =  0  (see 
appendix  C).  Therefore,  for  a  TEM  mode  to  propagate  in  the  sample,  the  constraint 
equations  for  the  air  gap  and  sample 


fe°)2+fe°)2  =  kl 

{K)2  +  (K)2  =  k] 
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- Fully  Filled  - 9  mil  Outer  Gap  -  No  Fix  - 9  mil  Inner  Gap  -  No  Fix 

- 9  mil  Outer  Gap  -  NIST  Correction  - 9  mil  Inner  Gap  -  NIST  Correction 


Frequency  [GHz]  Frequency  [GHz] 

(a)  (b) 


Figure  2.7:  Measured  permittivity  and  permeability  of  MagRAM  are  shown  as 

er  =  e'r  —  je"  and  fir  —  //fr  —  jfi",  repectively.  The  NIST  Correction  of  e'r  appears  to 
approach  the  actual  value  for  higher  frequencies,  but  does  not  perform  well  at  lower 
frequencies.  The  corrections  to  values  of  n'r  approach  the  actual  value,  but  are  not 
significant  enough.  Both  e"  and  /i"  are  biased  higher,  again  indicating  a  misreporting 
of  material  loss  properties. 
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- Fully  Filled  - 9  mil  Outer  Gap  -  No  Fix  - 9  mil  Inner  Gap  -  No  Fix 

- 9  mil  Outer  Gap  -  NIST  Correction  - 9  mil  Inner  Gap  -  NIST  Correction 


Figure  2.8:  In  figure  2.7(a)  the  corrected  e'r  of  MagRAM  appears  to  approach  the 
actual  value.  Zooming  in  at  2  GHz  shows  that,  while  realtive  close,  the  corrected 
values  for  both  inner  and  outer  gap  are  approximately  correct  and  nearly  follow  the 
slope  of  the  actual  value. 

are  reduced  to 


k°z  =  k0 
K  =  ks 


The  consequence  for  (2.2)  is 

e~jkoZ  =  e~jksZ  (2.3) 

must  hold  for  a  uniform  phase  front.  Of  course,  ko  =  cjyTo/Io  will  not  equal  ks  = 
/esfis  and  therefore  a  TEM  mode  cannot  exist. 

Furthermore,  two  assumptions  regarding  the  final  held  equations  are  made  to 
match  the  physical  expectation  of  the  system.  First,  the  higher  order  modes  excited 
by  the  air  gaps  in  Region  II  will  tend  to  match  the  symmetry  of  the  air  gaps.  Provided 
the  material  in  Region  II  is  axially  symmetric,  only  axially  symmetric  higher  order 
modes  are  expected. 
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Figure  2.9:  Considering  only  a  forward  propagating  TEM  excitation  and  ignoring 
any  reflections,  the  fields  in  Region  II  must  be  continuous  across  the  radial  material- 
air  gap  boundary.  Since  e~^k°z  7^  e_jfcsZ,  a  propagating  TEM  field  in  Region  II  is  not 
physically  possible. 


Second,  TE  modes  will  not  be  excited  within  Region  II.  The  fields  of  the  TEM 
excitation  are  (see  appendix  F  ) 


En  =  T 


(2.4) 


m  = 


T 


■1 TEM 


-,-jkz 


(2.5) 


where  T  is  a  constant  with  respect  to  frequency  and  the  material  parameters.  Axially 
symmetric  TE  fields  in  general  are 


Ep  =  Ez  =  H<j)  =  0  7^  0,  Hp  7^  0,  Hz  7^  0 

For  a  TE  held  to  be  excited,  the  incident  TEM  held  components  must  excite  orthogo¬ 
nal  helds  (i.e.  EjEIXI  must  excite  UjE).  Using  the  previous  assumption  of  an  isotropic 
media,  this  type  of  excitation  is  not  possible. 
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Axially  symmetric  TM  field  equations  within  the  sample  and  within  the  air  gap 
of  Region  II  must  independently  satisfy  the  wave  equation,  constraint  equation,  and 
applicable  boundary  conditions.  The  TM  field  equations  derived  in  appendix  F  are 
axially  symmetric  for  m  =  0  and  will  satisfy  these  conditions,  but  the  fields  within 
the  material  region  and  the  air  gap  region  will  be  unique  such  that  (considering  only 
forward  propagating  fields) 
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where  yn  =  jkz  for  mode  n,  A  and  B  are  constants  where  the  subscripts  indicate  the 
mode  and  region  where  the  constants  are  defined,  and  for  each  gap  scenario 
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Boundary  conditions  at  the  inner  and  outer  conductors  require  the  tangential  electric 
field  at  the  surface  to  be  zero 


Ez(p  =  a,  (j),  z)  =  Ez(p  =  a,  0,  z)  =  0 


which  are  applied  to  solve  for  constants  A  and  B 


o+  >0  {kpina) 
h0nJo(kplna) 


R+  Y0  ( kp2nb ) 
2,0n  J0  (kp2rib) 


Substituting  this  result  into  the  field  equations  produces 
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where 


kpl^nBiQn 
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kp2nlnB20n 
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The  hnal  boundary  condition  to  satisfy  requires  the  tangential  fields  to  be  con¬ 
tinuous  across  the  material  to  air  gap  interface  or 


H^R-)  =  H^(R+) 
EZ(R~ )  =  EZ(R  +) 


Applying  these  conditions  produces  two  separate  equations  for  the  constant  B20n 
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-°2,0n 
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D  + 
-°2,0n 
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kP2n  [Jo  ( kp2nR )  ho  ( kp2nb )  —  J0  ( kp2nb )  U0  (/cp2rii?)] 


which  are  equated  to  produce  the  characteristic  equation  for  %  , 


&pi„^TMin  [Jo  (kpinR)  ho  (kplna)  -  J0  (kplna)  h0  ( kplnR )] 

[Ji  ho  (fcp2»  -  Jo  (fcp2»  n  (fcp2„R)]  - 

kP2nZTM2ri  [Jo  ( kp2nR )  h"0  ( kp2nb )  -  J0  (/cP2n&)  h"0  (fcp2„-R)] 

[Ji  (kPinR)  ho  (kpina)  —  J0  (fcpi„a)  hi  (&pinR)]  =  0  (2.7) 


where 


^Pin  =  7n  +  and  &p2n  =  fn  +  cu2e2/r2 


Using  (2.7),  the  first  5  mode  wavenumbers  at  2  GHz  in  a  test  fixture  with  a  Rexolite® 
sample  and  100  mil  inner  and  outer  air  gap  are  calculated  and  for  MagRAM  with  a 
9  mil  inner  and  outer  air  gap  and  compared  to  wavenumbers  in  the  fully  filled  case. 
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The  results  are  reported  in  table  2.2  and  2.3.  It  is  immediately  evident  that  the 
propagating  TEM  field  corresponds  to  a  propagating,  axially  symmetric  TMoo  field 
in  both  materials. 

The  NIST  correction  factors  rely  on  the  similarity  of  the  transverse  field  com¬ 
ponents  of  the  TEM  and  TM0o  fields.  Recall  that  the  correction  factors  are  based  on 
the  material  modelled  as  a  lumped  capacitor  and  inductor  with  a  TEM  field.  The 
capacitance  and  inductance  per  unit  length  are  defined  as 
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respectively.  The  final  form  of  the  correction  factors  completes  the  integration  assum¬ 
ing  the  transverse  field  components  behave  in  a  -p  fashion,  the  expected  behavior  for 
a  TEM  field.  If  the  TM0o  field  components  within  the  material  and  air  gap  region  do 
not  follow  the  -  behavior,  the  accuracy  of  the  correction  factors  will  be  compromised. 

The  TM0o  field  components  for  a  Rexolite®  sample  with  an  100  mil  inner  and 
outer  air  gap  are  plotted  against  their  respective  TEM  field  components  in  figure  2.10. 
The  He/,  component  is  continuous  across  the  air  gap  to  material  boundary  as  required 
by  the  tangential  boundary  conditions  at  the  material  to  air  gap  interface.  The  Ez 
field  component  is  non-zero  as  expected  for  a  TMZ  field.  A  discontinuity  equal  to  the 
magnitude  of  the  permittivity  exists  in  Ep  for  both  gap  scenarios.  This  is  expected 
from  the  constitutive  relationship  for  a  simple  media  which  state 


p  ■  D  air 

_  P  _  -p  sample 

eo  E™ 

_  sample 

E;h 

=  er£psample 

It  is  this  discontinuity  that  (2.1)  is  accounting  for.  For  both  gap  scenarios,  the  fields 
within  both  the  air  and  material  exhibit  a  near  -p  behavior.  This  fact  is  exploited 
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by  the  NIST  correction  factor  and  justifies  the  excellent  corrected  e'r  values  shown  in 
figure  2.6(a).  Again,  correction  to  e",  shown  in  figure  2.6(b),  is  based  on  the  same 
TEM  assumption,  not  accounting  for  the  loss  of  power  due  to  higher  mode  excitation, 
and  therefore  is  not  expected  to  be  as  accurate  as  the  correction  to  e'r. 

The  TM0o  held  components  for  a  MagRAM  sample  with  an  9  mil  inner  and 
outer  air  gap  are  plotted  against  their  respective  TEM  field  components  in  figure  2.11. 
The  properties  of  the  fields  with  respect  to  an  Ep  discontinuity,  Ez  presence,  and  H $ 
continuity  at  the  material  to  air  gap  boundary  remain  the  same  as  the  Rexolite®  case. 
The  field  profiles  though  do  not  exhibit  a  -p  behavior.  This  violates  the  expected  -p 
behavior  and  mathematically  increases  the  complexity  of  the  integrations  in  (2.8). 
The  performance  of  the  correction  factors  is  reduced  as  shown  in  figure  2.7  (zoomed 
in  figure  2.8  for  the  real  permittivity).  The  imaginary  part,  for  both  permittivity  and 
permeability,  are  over  corrected  again  due  to  power  lost  to  higher  mode  excitation. 

2. 6  Summary 

Correction  factors  derived  by  modelling  a  sample  plus  axially  symmetric  air 
gap  as  a  lumped  capacitor  and  inductor  model  perform  well  when  the  assumptions 
regarding  the  model  are  met.  It  was  shown  that  for  a  low  loss  dielectric  material, 
the  assumption  of  a  TEM-like  field  within  the  sample  region  with  air  gap  is  valid 
and  produces  excellent  results  for  the  real  permittivity.  For  heavily  loaded,  magnetic 
samples  the  fields  within  the  material  region  with  air  gap  appear  less  TEM-like  and 
therefore  the  correction  to  the  real  permittivity  and  real  permeability  are  less  accurate. 
In  both  material  cases,  corrections  to  the  imaginary  part  of  the  permittivity  and 
permeability  are  hampered  by  not  accounting  for  power  coupling  into  excited  higher 
order  evanescent  modes. 
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Table  2.2:  Comparison  of  Wavenumbers  in  Rexolite®  when  inner  and  outer  air 

gaps  are  present.  For  7  the  frequency  is  2  GHz. 


Rexolite®  Filled  Coax 

100  mil  Outer  Air  Gap 

Mode 

7 

kp 

Mode 

7 

kp 

K 

TEM 

0.0+66.6j 

O.O+O.Oj 

=> 

TMqi 

0.0+59.5j 

29.9+O.Oj 

0.0+42.3  j 

TM01 

277.9 

285.7 

TMo2 

232.4 

241.8 

236.2 

tm02 

571.1 

575.0 

TMo3 

561.0 

564.9 

562.5 

tm03 

861.0 

863.5 

TMq4 

894.6 

897.1 

895.6 

TMq4 

1150.0 

1151.9 

tm05 

1171.3 

1173.2 

1172.0 

Rexolite®  Filled  Coax 

100  mil  Inner  Air  Gap 

Mode 

7 

kp 

Mode 

7 

kp 

K 

TEM 

0.0+66.6j 

O.O+O.Oj 

TM0o 

0.0+54.8j 

37.8+0.0.7 

0.0+35.4j 

TM01 

277.9 

285.7 

TMqi 

245.1 

253.9 

248.6 

tm02 

571.1 

575.0 

=> 

tm02 

569.7 

573.6 

571.2 

tm03 

861.0 

863.5 

TMo3 

899.3 

901.8 

900.3 

TM04 

1150.0 

1151.9 

TMq4 

1171.5 

1173.4 

1172.3 

Table  2.3:  Comparison  of  Wavenumbers  in  MagRAM  (er  =  21.03  —  z(17.71), 
pr  =  4.289  —  z(0.4633))  when  inner  and  outer  air  gaps  are  present.  For  7  the 
frequency  is  2  GHz. 


MagRAM  Filled  Coax 

9  mil  Outer  Air  Gap 

Mode 

7 

kp 

Mode 

7 

kp 

k°p 

TEM 

179. 3+419. 8j 

O.O+O.Oj 

TM00 

183.4+403.4* 

122.9-10.6* 

184.3+401.6* 

TM01 

224. 3+335. 7j 

285.7+O.Oj 

TM01 

133.9+214.3* 

364.0-128.0* 

135.8+211.3* 

tm02 

461.7+163.1.7 

575.0+0.0.7 

tm02 

347.6+197.9* 

475.3-13.7* 

349.5+196.9* 

tm03 

781. 6+96. 3j 

863. 5-0. Oj 

tm03 

663.5+122.6* 

754.5+8.0* 

664.7+122.3* 

TM04 

1089.8+69. lj 

1151. 9-0. Oj 

TMq4 

973.0+86.0* 

1040.9+8.0* 

973.9+85.9* 

MagRAM  Filled  Coax 

9  mil  Inner  Air  Gap 

Mode 

7 

kp 

Mode 

7 

kp 

K 

TEM 

179. 3+419. 8j 

O.O+O.Oj 

TM0o 

189.2+388.H 

171.4-10.8* 

190.1+386.3* 

TMqi 

224. 3+335. 7j 

285.7+O.Oj 

TMqi 

114.7+210.6* 

364.2-140.4* 

116.5+207.3* 

tm02 

461. 7+163. lj 

575.0+0.0.7 

tm02 

362.1+195.2* 

487.0-9.5* 

364.0+194.1* 

tm03 

781. 6+96. 3j 

863. 5-0. Oj 

tm03 

673.1+120.3* 

763.4+7.4* 

674.4+120.1* 

TMq4 

1089.8+69. lj 

1151. 9-0. Oj 

TMq4 

980.4+84.9* 

1047.9+7.6* 

981.3+84.8* 
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Figure  2.10:  Fields  at  2  GHz  for  a  Rexolite®  sample  with  a  100  mil  air  gap.  The 
Rexolite®  sample  properties  are  er  =  2.525  —  z(0.0)  and  \ir  —  1  —  z(0.0).  For  the 
coaxial  test  fixture,  the  inner  and  outer  conductor  radii  are  0.331  mil  and  0.76025  mil 
respectively.  The  fields  are  normalized  to  a  TEM  electric  field  of  K 


TEM 


TMno  Outer  Gap - TMQ0  Inner  Gap 


Figure  2.11:  Fields  at  2  GHz  for  a  MagRAM  sample  with  a  9  mil  air  gap.  The 
MagRAM  sample  properties  are  er  =  21.03  —  i(17. 71)  and  =  4.289  —  i(0.- 4633).  For 
the  coaxial  test  fixture,  the  inner  and  outer  conductor  radii  are  0.331  mil  and  0.76025 
mil  respectively. 
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III.  Full  Wave  Modal  Method 


Accurate  material  characterization  measurements  where  air  gaps  exist  between  the 
sample  and  test  fixture  conductors  require  consideration  of  the  higher  order  modes 
excited  by  the  air  gap.  The  full  wave  modal  method  describes  the  relationship  between 
the  fields  within  Regions  I,  II,  and  III  of  figure  3.1  as  a  function  of  the  frequency  and 
material  permittivity,  e,  and  permeability,  /i.  The  analysis  is  directed  at  determining 
a  theoretical  value  for  the  scattering  parameters,  S^y  and  given  the  specific  gap 
scenario  and  measured  gap  size.  The  calculated  theoretical  values  are  compared  to 
the  experimentally  measured  values1  as 


c  thy 
'-’li 


(u;,e,/r)  -  S 


exp 

11 


uy 


<  T 


c  thy 
°21 


(uye,^t) 


q  exp 
°21 


uy 


<  T 


(3.1) 


A  Newton  root  search  over  the  parameters  e  and  /i  is  executed  to  achieve  the  selected 
threshold,  T  (i.e.  (3.1)  is  evaluated  with  different  values  of  e  and  /i  until  both  equations 
are  true).  The  values  of  e  and  p  used  in  the  final  iteration  are  assumed  to  be  the 
actual  values  for  the  material  being  measured. 


3.1  General  Modal  Analysis 

Modal  analysis  is  a  method  of  solving  a  linear  system  for  unknown  parameters 
by  considering  all  the  potential  solutions  or  modes  of  the  system.  Applied  to  a  guided- 
wave  problem,  the  modes  are  the  electric  and  magnetic  field  particular  solutions  of 
the  Helmholtz  wave  equation.  A  primary  advantage  of  modal  analysis  is  the  clear 
physical  picture  maintained  through  the  development  that  not  only  calculates  the 
unknown  parameters,  but  readily  reveals  the  reason  for  the  solution.  The  three  steps 
of  modal  analysis  are 

1.  Describe  the  potential  modes  of  the  system  (field  expansion), 

1As  cited  previously,  the  reverse  measurements,  S22  and  S12  may  be  used  interchangeable  with 
the  forward  measurements,  Sn  and  S21 
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2.  apply  system  boundary  conditions  at  interfaces  to  ensure  uniqueness  of  the 
solution, 

3.  and  apply  a  selected  testing  operator  to  generate  a  linear  system  of  equations 
that  is  well-posed. 

The  resulting  system  of  equations  are  in  the  form 

Ax  =  b  (3-2) 

where  the  unknown  vector  x  is  readily  solved  for  via  linear  algebra  techniques.  Phys¬ 
ically,  the  modal  analysis  matrix  (MAM),  A,  for  the  guided  wave  system  of  this 
research  describes  the  coupling  between  held  modes.  The  forcing  vector  b  describes 
the  system  excitation  and  the  unknown  vector  x  describes  the  interaction  of  held 
modes  at  transverse  system  boundaries,  or  more  specifically  the  reflection  and  trans¬ 
mission  coefficients.  The  needed  scattering  parameters,  and  S^,  are  extracted 
from  the  vector  x. 

3.1.1  General  Definition  of  System  Modes.  An  infinite  number  of  modes 
are  possible  given  the  general  solution  of  helds  within  a  guided  wave  system  and 
application  of  the  appropriate  boundary  conditions.  The  total  transverse  electric  and 
magnetic  held  at  any  point  in  the  system  is  described  as  the  superposition  of  each 
potential  mode 

OO 

E  f  =  X>Je»eT7”*  (3.3) 

n=  1 

oo  „ 

±anKeTlnZ  ,  K  = 

i  £ n 

n=  1 

(3.4) 

where  the  subscript  t  represents  only  the  helds  transverse  to  the  direction  of  prop¬ 
agation  ( z ),  7 „  =  jkz  is  the  z  directed  wavenumber,  and  Zn  is  the  wave  impedance 
of  mode  n.  Coefficients  represent  the  complex  amplitudes  associated  with  each 
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mode.  The  coefficients  are  unique  and  constant  2  for  each  material  region  the  held  is 
defined  in  and  for  each  mode.  The  term  e^7"2  describes  the  field’s  propagation  along 
the  z  axis.3 

The  held  vectors  en  and  hn  are  mode  dependent  vectors  describing  the  electric 
and  magnetic  held  components  transverse  to  the  direction  of  propagation.  Disconti¬ 
nuities  in  the  coaxial  line  are  assumed  to  be  planar  and  perpendicular  to  the  primary 
axis  of  the  coaxial  line.  Therefore,  only  the  transverse  components  are  retained  in 
order  to  satisfy  the  tangential  boundary  conditions  at  transverse  discontinuities.  The 
complex  amplitude  and  held  vector  terms  may  differ  mathematically  for  each  n  (i.e. 
TEM  or  TM2  helds)  and  may  differ  over  the  cross  section  of  the  coaxial  line  (i.e. 
different  air  gap  sizes).  Therefore,  the  general  descriptions  of  (3.3)  and  (3.4)  are  used 
to  simplify  the  subsequent  development. 

The  sign  reversal  with  respect  to  direction  of  propagation  for  the  magnetic  held 
of  (3.4)  is  necessary  for  the  mathematical  description  to  match  the  physical  power 
how  in  the  system.  The  time  average  power  how  in  the  direction  of  propagation  is 
given  by  [5]: 

Pe  =  \{ExH*-n)  (3.5) 

where  n  is  the  held  propagation  direction.  Using  the  held  descriptions  in  (3.3)  and 
(3.4)  evaluated  at  z  =  0,  the  power  how  is 
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2The  materials  for  this  analysis  are  assumed  to  linear,  homogenous,  and  isotropic  (i.e.  simple 
media)  and  therefore  each  a„  is  a  constant. 

3By  convention  (exp5^),  the  term  eZlz  implies  a  field  propagating  in  the  ±£  direction 
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Forward  propagating  fields  [h  =  z )  exhibit  an  expected  positive  power  flow,  Pe,  along 
z.  Reverse  propagating  fields  (n  =  —z)  exhibit  a  non-physical  reverse  power  flow. 
Therefore,  the  sign  of  the  magnetic  field,  hn,  is  made  equivalent  to  the  sign  of  the 
direction  of  propagation  to  ensure  the  system  matches  the  physical  expectation  of  a 
positive  power  flow  in  the  direction  of  propagation. 

The  general  field  definitions  of  (3.3)  and  (3.4)  are  applied  using  figure  3.1  as 
a  guide  to  define  the  possible  held  modes  of  the  system.  Physically,  the  coaxial  line 
system  is  excited  from  Region  I  by  a  single,  z  travelling  TEM  held  produced  by  a 
network  analyzer.  Therefore,  (3.3)  and  (3.4)  in  Region  I  are 

OO  \ 

a^e-71*  +  a“ene+7™2 

^  <  0  (3.6) 

a^hie-712  -  £  a"hne+7"2 

n=  1  / 

where  the  frequency  of  excitation  is  assumed  below  the  cutoff  frequency  of  the  next 
higher  mode  and  therefore  only  a  single,  forward  propagating  mode  is  excited.  It 
is  also  assumed  that  the  NWA  is  perfectly  matched  to  the  system  so  there  are  no 
rehections  from  the  excitation  port  of  the  NWA  (see  hgure  2.1)4.  Therefore,  only  one 
forward  propagating  wave  (the  excitation)  will  exist,  yet  an  infinite  number  of  reverse 
propagating  helds  (rehections  from  the  Region  I  to  Region  II  boundary)  will  exist 
where  all  but  the  n  =  1  held  will  be  evanescent. 

The  configuration  of  helds  in  Region  II,  although  possessing  the  same  general 
form  (3.3)  and  (3.4),  may  be  significantly  different.  Therefore,  the  general  forms  are 
defined  uniquely  for  the  helds  in  Region  II  as 

OO  OO  x 

£  b+ene~^z  +  £  6~ene+7™2 

>  0  <z<L  (3.7) 

E  K K.e~lnZ  ~  E  K h„e+7"2 

n= 1  n= 1  / 


4In  practice,  the  network  analyzer  will  not  be  perfectly  matched  to  the  test  system.  A  perfectly 
matched  measurement  device  is  a  good  assumption  though  if  a  good  calibration  of  the  system  is 

accomplished  prior  to  measurement. 
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2  a 


2b 


V  -A/WV- 
*V\ AAA* 


-A/WV*  b1 


<V  -A/VW* 


-xMAA- 


a2’-"’an 


b2+,...,bn+  b2-,...,bn<y 


Figure  3.1:  The  air  gap  in  Region  II  presents  a  discontinuity  that  will  excite  higher- 
order  modes.  The  bandwidth  of  the  measurement  is  chosen  so  that  only  one  mode 
(subscript  1)  propagates  and  any  higher  order  modes  evanescence  (subscripts  2,...,n). 
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recognizing  that  the  material  and  makeup  of  Region  II  may  produce  different  held 
vectors  and  complex  amplitudes  then  Region  I. 

Region  III  is  filled  with  freespace  as  in  Region  I.  As  in  Region  I,  the  NWA  is 
assumed  to  operate  at  a  frequency  where  only  a  single  propagating  mode  is  supported 
by  the  coaxial  line  and  that  the  network  analyzer  is  perfectly  matched  to  the  system. 
Therefore,  the  fields  for  Region  III  are  given  by 

OO  -v 

c+ene_77l^“Z/) 

’S1  >  L  <  z  (3.8) 

£  c+hne~7ri(z_i) 

n=l  ) 

where  only  forward  propagating  fields  are  present  and  where  all  but  the  n  =  1  held 
will  be  evanescent.  It  follows  then  that  the  n  —  1  is  the  only  held  measured  by 
the  NWA5.  The  complex  amplitude,  c+,  is  referenced  to  the  boundary  at  z  =  L  by 
including  a  phase  shift  of  L , 


c+e7n£g-7n2  —  c+g-ln^-L) 

in  the  held  definition.  This  is  done  to  simplify  subsequent  analysis. 

3.1.2  Application  of  Boundary  Conditions.  The  transverse  components  of 
the  electric  and  magnetic  helds  are  continuous  across  the  transverse  interfaces  between 
Region  I  and  II,  and  Region  II  and  III  as  shown  in  hgure  3.1  [3].  Boundary  conditions 

5Certainly  evanescent  fields  carry  power  that  may  reach  the  NWA,  but  the  attenuation  is  signifi¬ 
cant  provided  the  NWA  detectors  are  reasonably  far  from  the  material  being  measured.  It  is  assumed 
that  any  evanescent  field  power  at  the  NWA  detectors  is  significantly  below  detection  capability  of 
the  NWA  and  therefore  insignificant. 
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at  z  =  0  and  z  =  L  due  to  the  presence  of  the  sample  in  the  coaxial  line  are 


hxE1(z  =  0  )  =  fixEn(z  =  0+) 
hxHI(z  =  0~)  =  nxH,J(z  =  0+)  +  Js(z  =  0) 
n  x  En(z  =  L~)  =  fix  E 111  (z  =  L+ ) 
h  x  Hn(z  =  L~)  =  hxHin(z  =  L+)  +  Js(z  =  L) 


Assuming  the  sample  has  no  impressed  sources  present  is  non-conducting  (Js  =  0) 
and  that  h  —  z ,  the  boundary  conditions  are  simplified  to 


Ete  =  (T)  =  E"(z  =  0+) 
H*  (z  =  0”)  =  UItI(z  =  0+ ) 
E  i\z  =  L~)  =  E  [n(z  =  L+) 
H  ['(z  =  L-)  =  H l"(z  =  L+) 


A  system  of  linear  equations  is  formed  by  substituting  the  field  equations  defined  in 
(3.6),  (3.7),  and  (3.8)  into  the  boundary  conditions  as 


OO 

O-t el  +  anen 

77.— 1 


OO  OO 

'y  1  bn  +  'y  ^  bn  en 

n—  1  77.— 1 


OO 

a+hi  -  y^a7hn 

77=1 


77=1  77=1 


^6+ene  7i  +  y^fcnene7i  =  J^c+e 

77=1  77=1  77=1 


^  6+hne  7L  -  ^  K  hne7i  =  c+  h 

77=1  77=1  77=1 


Recalling  that  the  scattering  parameters  are  the  ratio  of  incident  to  reflected  energy, 
the  system  is  scaled  by  Ap.  This  translates  the  constant  coefficients,  a,b,  and  c,  to 

ai 
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the  scattering  parameters  per  mode  at  each  transverse  interface  or 


.  -  b « 

_  + 

^ n  _i_ 

°1 

al 

t  -  bA 

Ln  _i_ 

ai 

a+ 

where  Rn  are  the  reflection  coefficients  at  z  =  0,  rn  are  the  ratio  of  fields  reflected 
at  z  =  L  and  propagated  to  z  =  0  over  the  incident  field,  tn  are  the  transmission 
coefficients  for  the  z  =  0  boundary,  and  Tn  are  the  transmission  coefficients  through 
the  material.  It  is  evident  that  the  terms  needed  for  the  comparisons  in  (3.1)  are 


Ri 


T\ 


x 


(3.9) 


By  substituting  the  scattering  parameters  per  mode  into  the  linear  system  and  trun¬ 
cating  the  summations  to  N  for  the  sake  of  computability,  the  linear  system 


N 

Y'.  Rnen 


n=l 

N 


Y  Rnhn 

n=  1 


N 

'Y. 

n=l 

N 

+  Y  f-nK. 

n=  1 
N 

tnene~lnL 

n= 1 

N 

£  tn 

n= 1 


N 

^  n^n 

n= 1 
N 

Yj  rnK 

n= 1 
N 

Y  rnene^L 

n= 1 

N 

Y  rnhne^L 

n= 1 


N 

Y  Tn^n 

n= 1 
N 

Y  TnK 

n= 1 


— ei 


hi 

(3.10) 

0 

0 


is  now  a  set  of  4  equations  with  4IV  unknowns.  Although  a  very  clear  physical  picture 
is  drawn  by  the  system,  it  is  just  as  clear  that  the  system  is  unsolvable  for  N  >  1. 
For  N  =  1,  the  system  is  well-posed,  but  does  not  account  for  higher  order  modes. 


3.1.3  Select  and  Apply  Testing  Operator.  A  testing  operator  is  applied  to 
the  system  of  equations  to  produce  a  well-posed  system  (i.e.  4N  equations  by  4N 
unknowns).  The  selection  of  the  testing  operator  should  be  physically  meaningful, 
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not  alter  the  scattering  parameter  definitions,  and  not  introduce  excessive  complexity 
into  the  system.  It  is  also  advantageous  computationally  if  the  set  of  testing  operators 
is  orthogonal  to  a  large  set  of  the  field  equations. 

The  Region  I  fields  meet  the  criteria  set  for  testing  operators.  Coupling  between 
two  field  modes  is  determined  by  multiplying  the  two  fields  together  and  integrating 
over  the  cross-section  of  the  coaxial  line  as  shown  in  section  2.2.  The  scattering 
parameters  are  constant  over  the  surface  of  integration  and  are  therefore  not  altered 
by  applying  the  testing  operator.  This  is  true  for  the  fields  in  all  three  regions,  but 
the  fields  within  region  I  and  III  (see  figure  3.1)  are  identical,  orthogonal  for  m  7^  n 
(as  shown  in  appendix  D),  and  unaltered  for  different  air  gap  widths  and  material 
properties.  Additionally,  the  fields  in  region  I  and  III  are  defined  by  a  single  expression 
valid  over  the  entire  cross-section,  whereas  the  region  II  fields  are  defined  by  unique 
expressions  in  the  material  and  air  gap.  Using  region  II  fields  would  require  two 
integrations  in  application  of  the  testing  operators,  whereas  using  region  I  fields  as 
the  testing  operator  results  in  a  single  integration  (per  mode,  per  operator)  when 
testing  in  region  I  and  III. 

The  testing  operator  is  applied  term-by-term  to  the  linear  system  of  equations 


as 


J  em  •  {}dS  ( 


m  =  1,  2, . . . ,  N)  ...  Operator  for  E-field  (3.11) 


cs 


m  =  l,2,...,N)  ...  Operator  for  H-field  (3.12) 


cs 
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where  the  integration  CS  is  over  the  cross  section  of  the  coaxial  line.  The  resulting 
system  of  linear  equations  is  well-posed  and  in  the  desired  form,  (3.1),  as 


(  A 
A, 


mn 

mn 


7  7 


—Br 

Cm, 


-Bn 

-C„ 


0 

0 


D  p  7 nL  TD  p'YnL 

C  p—'ynL  — plnL 
^  m.n.  °  m.n  c- 


0 

0 

— A„ 
Ar 


Zm'Z>n  / 


(  R  \ 
1  Ln 

tn 

rn 

\  Tn  ) 

\ 

zxz 


■Ami  ^ 

■Aml 
0 
0/ 


(3.13) 


where  rows  are  indicated  by  m,  columns  by  n,  and  the  submatrices  Amn,  Bmn,  and 
Cmn  are  defined  as 


A 

mn 

/  ® n  dS 

J 

(3.14) 

cs 

Bmn 

/  ■  en  dS 

J 

(3.15) 

cs 

C^mn 

=  f  hm  ■  hn  dS 

(3.16) 

cs 


The  relationship  between  e  and  h  eliminates  the  need  to  define  a  submatrix  for 
applying  a  magnetic  testing  operator,  h,  to  the  Region  I  and  III  magnetic  fields.  The 
electric  and  magnetic  transverse  mode  vectors,  en  and  hn  are  related  as 
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where  Zn  is  the  wave  impedance  for  mode  n.  Using  this  relationship,  the  product  and 
integration  of  the  testing  operator  and  fields  is  reduced  to 


hm  •  h ndS  = 


z  x  em  z  x  er 


dS 


cs 


cs 


7  7 

ZJmZJn 


(z  x  em)  ■  (z  x  en)  dS 


cs 


7  7 


{( %  '  5)(em  •  en)  ( Z  •  G n)(®m  ' 


C5 


Zm.Zri 


G777,  '  G ,, 


CS 


This  reduction  is  only  valid  in  Region  I  and  III  because  the  fields  in  these  regions 
are  defined  independently  of  the  material  or  air  gap  in  Region  II  and  the  fields  are 
defined  over  the  entire  surface  of  integration.  If  the  dimensions  of  the  coaxial  test 
fixture  remain  fixed,  the  fields  in  Region  I  and  III  will  vary  only  with  respect  to 
frequency.  This  reduction  reduces  the  number  of  unique  submatrices  required  in 
(3.13)  and  thereby  reduces  the  computational  load  in  forming  the  MAM. 

The  ill-posed  system  of  equations,  (3.10),  is  now  a  well-posed  system,  (3.13), 
that  generally  describes  the  inter-relationship  of  fields  in  the  three  coaxial  test  fixture 
regions.  The  system  is  in  the  form 


Ax  =  b  (3-17) 

where  the  vector  x  is  readily  solved  for  via  linear  algebra  techniques  and  contains  the 
scattering  parameters,  and  S^f  ,  necessary  to  complete  the  minimization  in  (3.1) 
and  determine  the  sample  permittivity,  e,  and  permeability,  ^t. 
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3.2  Application  of  Modal  Analysis 

The  modal  analysis  method  to  this  point  is  not  particular  to  a  coaxial  test 
fixture.  The  system  of  (3.13)  is  applicable  to  material  measurements  in  any  typical 
closed  waveguide  (i.e.  rectangular  or  cylindrical).  Although  this  research  considers 
an  air  gap  discontinuity,  this  fact  is  not  used  in  developing  the  system  of  (3.13).  The 
field  equations  within  Region  II  characterize  the  effects  of  the  air  gaps.  Therefore, 
to  specify  the  general  solution  to  a  particular  problem  and  waveguide  configuration 
requires  specifying  the  transverse  electric  and  magnetic  fields,  e„,  hn,  en,  and  hn,  in 
the  individual  regions  of  the  waveguide.  Following  this,  the  integrals  (3.14),  (3.15), 
and  (3.16)  must  also  be  evaluated  based  on  the  defined  transverse  mode  vectors  and 
air  gap  scenario  (inner  or  outer). 

3.2.1  Mode  Dependent  Transverse  Electric  and  Magnetic  Field  Vectors. 
Region  I  and  III  are  freespace-filled  coaxial  line.  The  fields  are  well-documented  [10] 
and  are  derived  in  appendix  F.  As  stated  earlier,  the  system  excitation  is  TEM  and 
the  axially  symmetric  air  gaps  in  Region  II  are  expected  to  excite  axially  symmetric, 
(m  —  0),  TM2  modes.  The  mode  dependent  transverse  field  vectors,  en  and  hn,  as 
defined  (3.6),  (3.7),  and  (3.8)  are  the  z-independent  terms  of  the  field  equations. 

Considering  only  forward  propagating  fields,  the  general  equations  in  appendix  F 
reduce  to 


EP  = 

= 


e~jkz 


_P 

1 


1 

P 


e~jkz 


(3.18) 

(3.19) 
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for  the  TEM  mode  and 


En  = 


E,  = 


H,  = 


TTjkp 

j2k 

7 Tk2 


Vi  ( kPnP )  e~lz 


j2k 
VT7  kp 


Pn  Vo  (kPnp)  e~lz 

“  1 


:ftk 


Hz  =  0 


Vi  (kPnp)  e~iz 

E(p  —  0 


Hp  =  0 


where 


Vm  ( ap )  =  Jm  ( ap )  Y0  (aa)  -  J0  (aa)  Ym  ( ap ) 


and  the  characteristic  equation  for  kPn  is 


Vo  (kPnb)  =  Jo  ( kPnb )  Y0  ( kPna )  -  J0  ( kPria )  F0  (kPnb)  =  0  (3.20) 

for  the  TMgn  modes.  Note  that  Vo  '  =  —V\  and  that  the  fields  above  are  normalized 
to  a  TEM  held  of  ke~jkz  using  the  constant  T,  defined  in  appendix  F  as 

Y  _ 

Tccope 

This  will  simplify  subsequent  calculations.  Based  on  these  held  definitions,  the  mode 
dependent  transverse  held  vectors  for  Regions  I  and  III  are 

(  1  - 

-  p  . . .  n  =  1 

P 

< 

'K^lnkp  / .  \  ~ 

■o,  Vd  (kPnp)  p  ...  n>  1 
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where 


h„  = 


1 

Zn 


7T7 nk 


Pn 


j2k 


Vl  (kPnP)  <t> 


n  =  1 


n  >  1 


7*  =  Kn-k 


k  =  ujy/JIe 


ry  _  7 n 

^ n  — 

JUJ6 

and  kPn  is  the  n  —  1  zero  of  (3.20).  For  n  —  1,  the  TEM  field,  kPl  =  0  and  therefore 
Z\  =  which  is  the  expected  wave  impedance  for  a  TEM  field. 

The  fields  in  Region  11  are  defined  in  section  2.5.  The  forward  propagating  field 
equations,  normalized  to  a  TEM  field  of  ^e-712  using  the  constant  T  just  as  in  the 
Region  1  and  111  fields,  are 
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(  nkPln% 

:ftk\ 
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^  1  (kplnp)  E\ne 


—InZ 


Wi  ( kPlnp )  Flne 
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“7  nZ 


Wi  ( kPlnp )  Flne 


—InZ 
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where 


kL  -  (^)2 


*L  - 


72 

in 


fn 


k\  = 

k2  — 

ll‘2e'2 

Zln  — 

In 

jve  i 

Z2n  — 

In 

jue  2 

wm  ( ap )  =  Jm  ( ap )  Y0  ( ab )  -  J0  (ab)  Ym  ( ap ) 


As  in  section  2.5,  the  boundary  conditions  at  the  material  to  air  gap  interface  requires 
that  the  tangential  fields,  Ez  and  H^,  are  continuous.  Equating  the  fields  at  p  =  R 
produces  two  continuity  factors  for  the  tangential  fields 


EZ(R  ~)  =  EZ(R  +) 


H^R  ~)  =  H^R  +) 


produces 

produces 


ln  ,n  VW 

p  _  p  kP2n  k  1  Zln  11  1  {kp2nR) 

ln  -  ln  kpin  k2  Z2n  V L  (fcpini?) 


and,  by  equating  the  continuity  factors,  the  characteristic  equation 


kPlnZ In  V0  (kpmR)  W1  (kP2nR)  ~  kP2nZ2n  K  (k^R)  W0  ( fcp2n  ^ )  =  0  (3-21) 

where  kPln  and  kP2n  are  the  nth  solution.  This  is  equivalent  to  (2.7),  but  is  reprinted 
here  in  simplified  form. 

Either  of  the  continuity  factors  can  be  inserted  into  the  field  equations  because 
the  characteristic  equation  ensures  they  are  equal  in  value.  Using  the  continuity 
factor  derived  from  is  instructive  though  in  that  the  ratio  of  the  impedances 
reduces  to  a  ratio  of  the  permittivities,  (A  Thus,  the  continuity  factor  readily 
shows  a  discontinuity  with  respect  to  the  permittivities.  This  is  expected  given  the 
constitutive  relationship  for  simple  media  (see  section  2.5). 
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The  continuity  factors  are  defined  above  to  substitute  for  Ein,  but  can  easily  be 
redefined  to  substitute  for  Fin.  The  choice  is  arbitrary.  For  this  research,  the  choice 
is  based  on  the  air  gap  scenario,  inner  or  outer.  Rather  than  define  simple  general 
equations  and  in  order  to  maintain  consistency,  the  continuity  factor  is  substituted 
into  field  equations  for  the  air  gap  in  Region  II  (i.e.  held  equations  defined  in  a  < 
p  <  R  for  inner  gap  or  R  <p<b  for  outer  gap). 

The  final  held  equations  for  an  inner  air  gap,  using  the  dimensions  defined  in 
figure  3.1,  are 


KkP2n%  £2  Wi  ( kP2nn ) 
j2k2  ei  Vi  (kPlnr i) 


Vi  ikpinp)  Flne  lnZ  ...  a  <  p  <  rx 


<  p  <  b 


'^kP2n'yn 

j2k2Z2n 


a  <  p  <  r  i 


H(tl  = 


( kP2np )  Flne~^ 


r  \  <  p  <  b 


a  <  p  <  r  i 


rx  <  p  <  b 
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Based  on  these  field  definitions,  the  mode  dependent  transverse  field  vectors  for  Region 
II  when  an  inner  air  gap  is  present  are 


e„  = 


^P2n7n  C2  lhj  {kP2nr i) 
.7  2 £>2  ci  hi 
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ri  <  p  <  b 


The  final  field  equations  for  an  outer  air  gap,  using  the  dimensions  defined  in 
figure  3.1,  are 
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Based  on  these  field  definitions,  the  mode  dependent  transverse  field  vectors  for  Region 
II  when  an  outer  air  gap  is  present  are 
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3.2.2  Mode  Coupling  Integrals  in  Coaxial  Line.  The  general  integrals  (3.14), 
(3.15),  and  (3.16)  can  now  be  defined  and  evaluated  based  on  the  mode  vectors  of  the 
previous  section.  The  solutions  are  elements  in  the  matrix  (3.13)  and  are  necessary 
to  solve  the  system  of  equations. 

The  surface  of  integration  in  each  case  is  the  cross-section  of  the  coaxial  line. 
Hence,  a  double  integral  with  respect  to  dp  and  pdcf)  are  expected  for  cylindrical  coor¬ 
dinates.  Axially  symmetric,  or  (j)  invariant,  field  excitation  is  a  primary  assumption  in 
the  preceding  development  and  therefore  the  integration  over  dcj)  reduces  to  a  constant 
27t  factor.  For  example,  (3.14)  becomes 
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2tt  em-  enpdp 


The  factor  of  2n  is  present  in  all  the  integrations  necessary  to  populate  (3.13)  and  is 
therefore  normalized  out  of  the  following  development. 
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Coupling  between  fields  in  the  freespace  region  (Region  I  and  III)  are  determined 
from  (3.14).  Knowing  that  application  of  the  testing  operator  is  identical  in  form  to 
testing  for  orthogonality  and  that  the  transverse  fields  in  this  region  are  orthogonal 
(see  appendix  D),  the  solution  to  (3.14)  is  0  for  m  ^  n.  The  self-coupling,  m  =  n, 
solutions  to  (3.14)  are  in  two  forms.  The  first  for  the  TEM,  m  =  n  =  1,  field 
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and  the  second  for  the  TM,  m  =  n  >  1,  fields 
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The  integration  is  accomplished  using  the  identity  (E.l)  and  produces 
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which  is  greatly  simplified  to 
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by  recognizing  the  following  equalities 


MkPna) 


Vo(kPnb ) 


Vl  (kPna) 


0  ...  trivial 

0  ...  TM  field  characteristic  equation  (3.20) 

-  ...  Wronskian  relationship  [1] 


The  remaining  mode  coupling  integrations,  (3.15)  and  (3.16),  require  lengthy 
proofs  to  simplify.  General  observations  and  the  final  solutions  are  provided  here,  but 
the  proofs  are  relegated  to  appendix  G. 

Inner  and  outer  air  gap  scenarios  are  shown  in  figure  3.1.  The  air  gap  disconti¬ 
nuities  are  invariant  in  (f)  and  therefore  integration  with  respect  to  d(j)  is  a  constant, 
27T,  as  previously  noted  and  is  normalized  out.  The  fields  in  the  air  gap  and  material 
regions  are  not  equivalent  and  therefore  the  discontinuity  in  p  requires  the  general 
integrations  of  (3.15)  and  (3.16)  be  split  to 
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(3.27) 

(3.28) 


where  R  —  r\  for  an  inner  air  gap  and  R  =  r 2  for  an  outer  air  gap.  Furthermore,  there 
are  two  distinct  forms  of  the  testing  operators  em  and  hm,  namely  TEM  (m  =  1)  and 
TM  (m  >  1),  and  two  distinct  forms  of  the  field  equations  en  and  h„,  one  set  of 
equations  given  an  inner  air  gap  and  a  different  set  given  an  outer  air  gap. 

Close  comparison  of  the  inner  and  outer  air  gap  transverse  mode  vectors  in 
Region  II,  (3.22)  to  (3.24)  and  (3.23)  to  (3.24),  reveals  that  variations  with  respect  to 
p  are  identical  for  each  set  of  equations.  Only  the  constants  change  depending  on  an 
inner  or  outer  air  gap  scenario.  This  is  true  for  both  the  electric  and  magnetic  fields 
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because 


where  Zn  is  independent  of  p.  Thus,  rather  than  the  four  integrations  in  (3.27)  and 
(3.28)  being  carried  out  for  inner  and  outer  gap  scenarios,  and  for  a  TEM  and  TM 
testing  operator  (a  total  of  16  unique  integrations),  only  4  unique  integrations  are 
necessary  with  respect  to  the  p  varying  terms  and  then  modified  by  the  appropriate 
constants.  This  not  only  simplifies  the  development,  but  significantly  reduces  the 
computation  time. 

These  insights  permit  simplifying  the  general  proofs  of  appendix  G  to  compu¬ 
tationally  efficient,  air  gap  scenario  (inner  or  outer)  dependent  solutions.  The  result 
for  an  outer  air  gap  scenario  is 
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and  for  m  ^  1 
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The  result  for  an  inner  air  gap  scenario  is 
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and  for  m  ^  1 
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The  integrations  necessary  to  fill  the  matrices  of  (3.13)  for  both  an  inner  and 
outer  air  gap  scenario  are  complex,  but  computationally  possible.  Care  must  be 
taken  in  determining  the  values  of  kpm  from  (3.20),  and  kp \n  and  kp2n  from  (3.21). 
Additionally,  if  the  values  of  k2rn  —  kp2n  or  kpm  —  kpln  approach  0  the  form  of  (3.36), 
(3.37),  (3.31),  and  (3.32)  will  be  unstable. 


3. 3  Summary 

A  method  for  calculating  theoretical  scattering  parameters  was  developed  for 
material  characterization  measurements  where  the  material  does  not  completely  fill 
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the  cross  sectional  area  of  a  coaxial  line  test  fixture.  Expressions  for  filling  the  elements 
of  the  MAM  are  given  for  both  an  inner  and  outer  air  gap.  After  filling  the  MAM 
and  forcing  vector,  b,  solving  for  the  reflection  and  transmission  coefficient  vector 
in  (3.13)  is  simply  a  matter  of  linear  algebra.  The  needed  terms  and  are 
easily  extracted  from  the  solution  vector,  x,  and  compared  to  the  measured  scattering 
parameters.  If  the  difference  in  theoretical  and  measured  scattering  parameters  is 
below  a  user  provided  tolerance,  the  permittivity  and  permeability  used  in  the  MAM 
calculations  are  representative  of  the  material.  If  not,  a  Newton  root  search  iterates 
the  values  of  permittivity  and  permeability  until  the  desired  tolerance  is  met. 
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IV.  Full  Wave  Modal  Method  Results 


The  Full  Wave  Modal  Method  was  developed  in  the  previous  chapter  for  the  cases 
of  a  sample  with  inner  or  outer  air  gap.  The  objective  of  the  analysis  is  to  calcu¬ 
late  theoretical  scattering  parameters  S^y  and  which  are  then  compared  against 
experimentally  measured  scattering  parameters.  The  values  of  permittivity,  e,  and 
permeability,  /i,  for  a  given  frequency,  u ;,  are  iterated  until  the  difference  is  below  a 
specified  tolerance,  T,  or 
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and  the  final  values  of  e  and  p  are  that  of  the  measured  sample. 

The  algorithm  uses  a  Newton’s  root  search  to  accomplish  the  minimization  in 
4.1.  For  a  non-magnetic  material,  only  a  ID  complex  root  search  is  needed  (/ir  =  1), 
whereas  a  2D  complex  root  search  is  used  for  magnetic  material.  In  either  case,  the 
algorithm  requires  an  initial  guess  that  is  updated  with  each  iteration  until  the  final 
solution  is  reached. 


4-1  Test  Setup 

The  method  described  in  chapter  111  is  implemented  in  Matlab®.  Measurements 
of  the  Rexolite®  and  MagRAM  materials  shown  in  chapter  11  are  evaluated  using 
the  modal  method  and  the  results  reported  herein.  Although  multiple  samples  of 
Rexolite®  and  MagRAM  were  used,  all  of  the  samples  are  assumed  to  be  homogeneous. 
The  material  dimensions  and  electromagnetic  properties  are  described  in  section  2.3. 

The  coaxial  test  fixture  is  shown  assembled  in  figure  4.2.  The  coaxial  line  inner 
conductor  radius  is  0.331  inches  and  outer  conductor  radius  is  0.76025  inches,  as 
described  in  chapter  11.  The  test  fixture  center  conductor  is  not  held  together  with 
hardware.  The  sample  holder  center  conductor  (detail  shown  in  figure  1.1  and  4.1) 
is  flanged  on  both  ends.  The  center  conductor  in  the  sample  holder  is  tension  fit 
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Figure  4.1:  Coaxial  test  fixture  opened  with  center  conductor  inserted. 

into  mating  center  conductor  pieces  that  run  the  length  of  the  test  fixture  shown  in 
4.2.  The  flanges  are  soft,  but  are  adjusted  when  removed  to  ensure  a  good  electrical 
contact  is  made  between  the  two  surfaces  when  reconnected. 

The  coaxial  test  fixture  is  connected  to  a  Network  Analyzer  for  measuring  the 
scattering  parameters  from  50  MHz  to  3.05  GHz  as  described  in  chapter  II.  As 
mentioned  previously,  the  system,  NWA  and  test  fixture  (including  the  sample  holder), 


Figure  4.2:  Assembled  coaxial  test  fixture. 
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are  calibrated  via  the  TRL  method  prior  to  measurement  and  care  is  taken  to  limit 
cable  movement. 


4-2  Error  Analysis 

Sources  of  potential  error  include,  but  are  not  limited  to,  uncertainty  in  sample 
length  measurements  and  radial  thickness  of  sample.  All  measurements  were  made 
with  a  calibrated  set  of  calipers  with  precision  of  0.5  mils.  Repeated  measurements  by 
a  second  person,  showed  general  accuracy  of  ±  1.0  mil  from  the  originally  measured 
values1.  Both  the  sample  length  and  radial  thickness  dimensions  were  measured  using 
the  same  method  and  hence  the  observed  uncertainty  in  measuring  is  assumed  to  be 
equally  applicable  to  each  source. 

An  error  analysis  for  the  real  and  imaginary  permittivity  and  permeability  are 
accomplished  as  in  [7].  The  permittivity  and  permeability  as  a  function  of  sample 
length,  L ,  and  sample  radius  thickness,  R.  are  expanded  in  a  Taylor  Series  around 
the  original  measurements  L0  and  i?0  plus  measurement  uncertainty  A L  and  A R, 
respectively.  The  modal  method  is  not  differentiable,  as  expected  by  the  Taylor 
Series,  but  if  higher-order  terms  are  assumed  negligible  and  the  errors  assumed  to 
be  relatively  small,  the  exact  error  is  approximated  by  a  sum  of  partial  derivatives, 
shown  here  for  the  real  permittivity,  as 
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lrThe  radial  thickness  of  several,  random  samples  were  measured  a  second  time  by  a  second, 
less-experienced  individual.  The  first  measurements,  made  by  an  individual  experienced  with  the 
equipment,  is  assumed  to  be  the  actual  value. 
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which  are,  in  turn,  approximated  using  finite  difference  derivatives.  The  result,  again 
shown  for  the  real  permittivity  is  a  simple  form 
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that  individually  accounts  for  sources  of  error.  A  worst  case  error  estimate  is  achieved 
by  a  composite  of  the  error  quantities  as  a  triangle  inequality 


\^e'approx\  -  |Ae'|  -  |Ae  'L  +  Ae  ^|  <  |Ae  'L\  +  |Ae^| 

where  the  magnitude  of  the  individual  terms  is  used  to  prevent  error  cancellation 
when  the  terms  are  opposite  in  sign.  Assuming  the  errors  AL  and  A R  are  relatively 
small,  the  error  is  assumed  to  be  symmetric  around  the  actual  value  (i.e.  +2  mil 
measurement  error  will  produce  the  same  results  as  -2  mil  measurement  error).  Error 
bars  on  graphs  are  a  result  of  analyzing  uncertainty  in  sample  length  and  radial 
thickness  of  sample  measurements  of  ±2  mil  measurement  to  produce  a  conservative 
estimate  of  the  error. 


4-3  Modal  Method  Performance 

The  same  NWA  measurements  used  to  produce  figures  2.4  and  2.5  for  Rexolite® 
with  100  mil  inner  and  outer  air  gap  and  MagRAM  with  9  mil  inner  and  outer  air  gap 
respectively,  are  now  used  here  in  evaluating  the  performance  of  the  modal  method. 

4-3.1  Modal  method  performance  for  Rexolite®.  The  modal  method  per¬ 
formed  very  well  for  the  Rexolite®  data  considering  10  modes.  The  permittivity,  e'r 
and  e",  is  extracted  with  a  high  degree  of  accuracy  as  shown  in  figure  4.3.  As  a 
low-loss,  low  dielectric  constant  material,  the  wavenumbers  for  Rexolite®  are  easily 
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predicted,  producing  a  single  propagating  mode  and  9  highly  evanescent  modes  (see 
table  2.2). 

The  results  for  e"  in  Rexolite®  are  as  expected  and  show  very  little  error.  Al¬ 
though  the  e"  are  negative  at  lower  frequencies,  the  magnitude  of  the  values  (-.015 
being  the  peak)  suggests  noise  in  the  system  or  other  unaccounted  for  anomaly,  rather 
than  calculation  error.  There  are,  however,  two  noticeable  phenomenon  in  the  e'r  re¬ 
sults,  namely  inner  air  gap  error  and  data  oscillations. 

The  inner  air  gap  calculation  for  e'r  is  not  as  accurate  as  the  outer  air  gap  when 
considering  the  same  number  of  modes.  For  a  coaxial  line  filled  with  freespace  and 
considering  the  electric  field  component  of  a  TEM  wave,  32%  of  the  potential  is  in 
the  space  between  the  inner  conductor  and  a  point  100  mil  out  radially.  Only  17% 
of  the  potential  is  contained  in  the  space  between  the  outer  conductor  and  a  point 
100  mil  in  radially.  Applied  to  the  air  gap  scenario,  the  sample  material  perturbs 
less  of  the  electric  field  strength  when  an  inner  air  gap  is  present  versus  an  outer  air 
gap.  Therefore,  analysis  where  an  inner  gap  is  present  is  more  susceptible  to  error, 
including  that  from  uncertainty  in  the  radial  width  measurement  (as  shown  in  figure 
4.4). 

Oscillations  in  the  results  are  attributed  to  effects  within  the  test  system  that 
were  not  accounted  for  in  the  calibration.  A  TRL  calibration  of  the  NWA  and  coax¬ 
ial  test  fixture  requires  breaking  the  inner  conductor  two  times.  The  measurement 
process  requires  the  coaxial  line  be  broke  again  to  insert  the  test  sample.  The  act 
of  breaking  and  reassembling  the  coaxial  line,  even  during  calibration,  will  effect  the 
quality  of  the  calibration  and  ultimately  the  measurement.  The  oscillations  present 
in  figure  4.3  are  not  significant  enough  to  cause  concern. 

The  impact  of  sample  length  measurement  on  the  calculation  of  e"  (figure  4.4(b)) 
is  nearly  imperceptible,  but  worth  noting.  Assuming  a  PEC  coaxial  line,  loss  within 
the  sample  material  will  be  different  based  on  the  amount  of  material  the  fields  travel 
through.  Therefore,  uncertainty  in  the  sample  length  measurement  will  have  a  greater 
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impact  on  the  the  calculation  of  e"  versus  uncertainty  in  the  radial  width  measurement 
of  the  sample.  The  small  magnitude  of  the  error  in  e"  is  a  reflection  of  uncertainty, 
±2  mil,  relative  to  the  wavelength,  3.934  in.  at  3  GHz. 

A  primary  benefit  of  the  modal  method  is  the  convergence  of  the  result  to  the 
material  properties.  As  shown  in  figure  4.5,  using  additional  modes  in  the  analysis 
does  increase  the  accuracy  of  the  result  in  an  expected  manner,  but  the  accuracy  gain 
is  insignificant  beyond  10  modes.  Additionally,  increasing  the  number  of  modes  in  the 
analysis  only  compensates  for  higher  order  mode  excitation,  not  for  other  errors  in 
the  system.  Radial  width  or  sample  length  measurement  uncertainty  error  does  not 
decrease  with  an  increase  in  the  number  of  modes  in  the  analysis.  In  figure  4.5,  the 
number  of  modes  is  increased  to  25.  Error  due  to  the  radial  width  and  sample  length 
measurement  uncertainty  converges  at  the  rate  the  calculated  results  do  proving  that 
the  introduction  of  additional  modes  only  compensates  for  error  due  to  the  excitation 
of  those  higher  order  modes. 

Modal  method  performance  for  different  gap  sizes  was  tested  as  well.  In  total, 
24  Rexolite®  samples  were  machined  to  produce  outer  and  inner  gaps  of  1,  2.5,  3.75, 
5.25,  20,  40,  60,  80,  100,  200,  300,  and  400  mils.  As  seen  in  figure  4.7,  the  modal 
method  is  able  to  accurately  calculate  the  permittivity  for  each  outer  air  gap  size 
(data  shown  is  for  measurements  at  2  GHz).  The  modal  method  applied  to  inner 
air  gaps  (not  shown),  produced  a  similar  result  up  to  300  mil.  Not  surprisingly,  the 
400  mil  inner  air  gap  results  are  poor.  For  a  400  mil  inner  air  gap,  only  29.25  mil  of 
Rexolite®  lines  the  inside  of  the  outer  conductor  and  only  5%  of  the  electric  potential 
impinges  on  the  material. 

The  lumped  circuit  correction,  developed  in  chapter  11,  is  compared  to  the  modal 
method  for  varying  gaps  in  figure  4.7.  The  lumped  circuit  model  provides  good  results 
up  to  a  60  mil  outer  air  gap  and  its  behavior  for  greater  gap  sizes  is  similar  to  the 
modal  method  although  different  in  magnitude.  The  lumped  circuit  correction  is 
applied  to  the  NRW  algorithm  output.  The  NRW  algorithm  uses  the  same  measured 
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data  that  the  modal  method  does,  but  assumes  only  single  mode  propagation.  For 
that  reason,  the  NRW  algorithm  and  lumped  circuit  model  behavior  will  compare  well 
to  the  modal  method  when  a  finite  number  of  modes  are  used  and  provides  a  good 
confidence  metric  for  the  modal  method. 

A  critical  assumption  of  the  above  theory  is  an  axially  symmetric  air  gap.  Fail¬ 
ing  this  assumption  will  not  produce  the  expected  results,  even  for  small  amounts 
of  conductor  or  material  eccentricity.  In  the  case  of  a  fully-filled  coaxial  line  mea¬ 
surement,  the  optimal  situation,  the  center  conductor  is  physically  supported  by  the 
sample  material.  When  an  air  gap  is  introduced,  proper  support  must  be  added  to 
the  center  conductor  to  ensure  the  weight  of  the  center  conductor  (for  an  inner  air 
gap)  or  the  weight  of  the  center  conductor  plus  sample  (for  an  outer  gap)  do  not  pull 
the  center  conductor  off  the  central  axis  of  the  coaxial  line. 

A  test  system  is  easily  evaluated  for  an  eccentric  center  conductor  by  measuring 
a  known  sample  with  known  air  gap  and  comparing  the  results  as  the  number  of  modes 
is  increased.  An  example  is  shown  in  figure  4.6,  where  a  sample  of  Rexolite®  only 
29.25  mil  thick  (producing  a  400  mil  outer  air  gap)  is  analyzed  using  20  modes.  For 
the  eccentric  case,  the  center  conductor  was  measured  as  6  mil  off  the  center  axis  of 
the  coaxial  line.  The  eccentricity  of  the  center  conductor  limits  the  convergence  of  the 
permittivity  (both  real  and  imaginary)  to  a  value  different  than  the  known  value  for 
the  sample.  Additionally,  the  results  for  e”  incorrectly  indicate  the  material  is  lossy, 
particularly  at  lower  frequencies.  Supports  for  the  center  conductor  were  than  added 
and  the  measurement  retaken.  The  modal  method  with  20  modes  was  used  again  to 
calculate  the  permittivity  and  the  results,  also  shown  in  figure  4.6,  match  the  expected 
value  very  well.  Additionally,  e"  converges  to  the  expected  value  over  frequency.  The 
low  frequency  values  of  e"  are  non-physical  and  noisy.  The  finite  conductivity  of  the 
coaxial  line  adversely  impacts  the  measurements  at  low  frequencies,  but  the  noise  in 
the  calculation  is  a  good  indication  (for  this  thin  of  a  sample)  that  the  major  source 
of  error,  higher  order  mode  excitation,  is  accounted  for. 
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Fully  Filled - 100  mil  Outer  gap  (10  Modes) - 100  mil  Inner  gap  (10  modes) 


(a)  (b) 

Figure  4.3:  Measured  permittivity  is  shown  as  er  =  e'r  —  je”.  Error  bars  indicate 

potential  error  in  radial  width  (±2  mil)  and  length  (±2  mil)  measurements.  Using 
the  modal  method  with  10  modes  considered  produces  a  very  accurate  correction  for 
the  real  (a)  and  imaginary  (b)  part  of  the  permittivity.  Oscillations  in  the  results  are 
on  the  order  of  .1  for  the  real  part  and  .01  for  the  imaginary  part.  The  oscillations 
are  attributed  to  system  imperfections  not  accounted  for  in  the  calibration. 

4-3.2  Modal  method  performance  for  MagRAM.  The  modal  method  per¬ 
formed  reasonably  well  for  the  MagRAM  data  with  a  9  mil  inner  and  outer  air  gap 
considering  10  modes.  The  outer  air  gap  result,  shown  in  figure  4.8,  is  fairly  accu¬ 
rate  and  follows  the  behavior  of  the  fully  filled  measurements  for  all  parameters  with 
some  error  still  present  in  the  e”  result.  The  inner  air  gap  data  does  not  follow  the 
behavior  of  the  fully  Elled  data  very  well  at  low  frequencies  for  permittivity  and  high 
frequencies  for  permeability. 

The  error  evident  in  the  MagRAM  permittivity  results  indicates  one  of  four 
things.  First,  and  most  likely,  eccentricity  of  the  center  conductor  is  again  impacting 
measurements.  Although  only  a  small  gap  exists,  the  material  is  heavily  loaded.  The 
results,  particularly  e ",  indicate  that  the  fields  in  the  system  are  not  exclusively  axially 
symmetric.  This  conjecture  requires  additional  measurements  to  prove. 

Numerical  instability  of  the  Bessel  functions  with  complex  arguments  is  consid¬ 
ered  the  second  most  likely  cause  of  error.  Bessel  functions  are  a  fundamental  piece 
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2  mil  radial  width  measurement  error 


2  mil  sample  length  measurement  error 


Frequency  [GHz] 


Frequency  [GHz] 


(a) 


(b) 


Figure  4.4:  Permittivity  of  Rexolite®  with  100  mil  inner  air  gap  using  10 

mode  analysis  is  shown  for  measurement  error  in  sample  radial  width  (air  gap  is 
100  mil  ±2  mil)  and  length  (sample  length  is  1  in  ±2  mil).  Measurement  of  e'r  is 
more  impacted  by  error  in  the  radial  width  measurement  because  of  the  high  field 
concentration  at  the  center  conductor.  Measurement  of  e"  is  impacted  more  by  the 
sample  length  measurement  error  as  a  homogeneous  sample  will  produce  more  or  less 
loss  in  the  system  as  the  length  is  decreased  or  increased  respectively.  Error  in  the 
real  (a)  and  imaginary  (b)  part  of  the  calculated  permittivity  is  inconsequential. 


Fully  Filled - 100  mil  Outer  gap  - 100  mil  Inner  gap 


Figure  4.5:  Measured  permittivity  at  2  GHz  shown  as  er  =  e'r  —  je".  Error  bars 

indicate  potential  error  in  radial  width  (±2  mil)  and  sample  length  (±2  mil)  measure¬ 
ments.  Increasing  the  number  of  modes  does  not  impact  the  error  from  radial  width 
and  length  measurments.  For  low  loss,  low  dielectric  constant  Rexolite®,  the  results 
converge  with  10  modes. 
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- TRUTH - 20  Modes  - 20  Modes  -  Eccentric  Center  Conductor 


(a)  (b) 

Figure  4.6:  Calculated  permittivity  for  a  Rexolite®  sample  29.25  mil  thick  (400  mil 
outer  air  gap)  using  20  modes  and  shown  as  er  =  e'r  —  je".  An  eccentric  center 
conductor  will  excite  non-axially  symmetric  higher  order  modes  that  are  not  accounted 
for  in  this  development.  Extra  support  was  added  to  ensure  the  center  conductor 
remained  along  the  center  axis  of  the  coaxial  line  so  that  only  axially  symmetric 
modes  were  excited. 


- Fully  Filled  -®-  20  modes  — —  Lumped  Circuit 


Figure  4.7:  Calculated  permittivity  for  a  Rexolite®  shown  as  er  =  e'r  —  je".  Mea¬ 
surements  are  at  2  GHz.  The  outer  air  gap  is  varied  (by  using  multiple  samples)  from 
1  mil  to  400  mil.  The  lumped  circuit  is  shown  for  comparison. 
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to  the  theory  presented  here.  Wavenumbers  for  higher  order  fields  in  the  Rexolite® 
samples  were  exclusively  real  due  to  the  measurement  frequency  band.  The  lossy 
nature  of  the  MagRAM  material  insures  the  higher  order  wavenumbers  are  complex 
and,  as  mentioned  previously,  the  real  part  grows  as  the  imaginary  part  shrinks  signif¬ 
icantly.  The  Bessel  functions  must  be  very  accurate.  As  seen  in  figure  4.15,  the  wave 
numbers  calculated  for  the  MagRAM  material  are  non-physical  beyond  25  modes.  A 
more  stable  Bessel  calculation  is  expected  to  address  this  potential  source  of  error. 
The  effects  are  discussed  further  in  section  4.4.2. 

Thirdly,  the  coaxial  is  not  a  perfect  conductor,  but  is  assumed  to  be  such.  If 
the  actual  conductivity  of  the  conductor  is  used  in  the  development,  the  boundary 
conditions,  and  therefore  the  fields,  are  different  than  those  used  in  this  research.  The 
effect  will  be  particularly  pronounced  at  lower  frequencies.  Exploring  these  effects  is 
worthwhile,  but  outside  the  scope  of  this  research. 

The  fourth  most  likely  cause  of  the  error  in  the  permittivity  is  inhomogeneity 
of  the  MagRAM  sample.  Care  was  taken  in  cutting  samples  from  MagRAM  sheets, 
but  it  is  not  unlikely  that,  with  such  a  heavily  loaded  material,  the  samples  may  be 
of  different  dielectric  composition. 

The  accurate  and  stable  measurement  of  the  permeability  is  particularly  no¬ 
tably.  An  outer  gap  scenario  is  the  most  accurate,  as  mentioned  previously,  but  both 
air  gap  scenarios  shown  in  figure  4.8  provide  good  results  for  permeability  over  the 
measurement  frequency.  Figure  4.10  indicates  the  permeability  results  are  also  very 
stable  with  respect  to  the  number  of  modes  included  in  the  analysis.  The  magnetic 
field,  Hip,  is  continuous  across  the  material  air  gap  boundary  as  shown  in  figure  2.11. 
This  is  required  by  tangential  boundary  conditions.  As  such,  the  H $  held  and  there¬ 
fore  the  results  are  not  as  sensitive  to  modest  variations  in  the  permeability. 

Error  caused  by  uncertainty  in  the  sample  length  measurement  is  largely  over¬ 
shadowed  by  error  from  uncertainty  in  the  radial  width  measurement  in  figure  4.9. 
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The  effect  is  identical  to  that  described  for  the  Rexolite®  material,  but  is  significantly 
enhanced  due  to  the  heavy  dielectric  loading  in  the  MagRAM  sample. 

The  permittivity  results  converge  fairly  rapidly  although  not  as  quickly  as  the 
Rexolite®  results,  with  the  largest  gains  within  the  first  10  modes  (figure  4.10).  The 
boundary  conditions  require  a  discontinuous  Ep  field  across  the  material  to  air  gap 
boundary  and  the  characteristic  equation,  (3.21),  is  largely  dependent  on  this  discon¬ 
tinuity.  Viewed  from  a  field  perspective,  the  higher  order  fields  that  are  excited  do 
not  decay  as  rapidly  in  the  heavily  loaded  MagRAM  as  in  the  Rexolite®  samples  (ev¬ 
ident  by  comparison  of  the  wavenumbers  in  table  2.2  and  2.3).  Hence,  the  inclusion 
of  higher  order  fields  is  important  to  accurately  determining  the  permittivity. 

Modal  method  performance  for  different  gap  sizes  was  tested  as  well.  In  total, 
8  MagRAM  samples  were  machined  to  produce  outer  and  inner  gaps  of  3,  6,  9,  and 
12  mils.  The  modal  method  is  able  to  accurately  calculate  the  permittivity  for  each 
outer  air  gap  size,  shown  in  figure  4.11.  The  modal  method  applied  to  inner  air  gaps 
(not  shown),  produced  similar  results  for  all  the  gap  sizes. 

The  lumped  circuit  correction  again  is  a  good  confidence  metric  for  the  modal 
method.  The  lump  circuit  correction  results,  compared  in  figure  4.11,  have  the  same 
behavior  and  similar  values. 

4-4  Root  Search  Algorithm  Performance 

This  research  uses  both  Newton  and  Muller  root  searches.  The  Newton  root 
search  essentially  follows  linear  paths  looking  for  the  smallest  value  on  a  line  (ID) 
or  surface  (2D).  The  Newton  root  search  only  requires  a  single  initial  guess  for  each 
variable.  A  Muller  root  search  is  similar  to  the  Newton  method,  but  uses  inverse 
quadratic  interpolation  rather  than  a  linear  interpolation  [18] .  The  Muller  root  search 
was  used  when  the  roots  of  an  equation  might  be  in  close  proximity  and  therefore 
small,  non-linear  interpolations  were  necessary  to  ensure  the  appropriate  root  was 
found. 
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- Fully  Filled - 9  mil  Outer  gap  (10  Modes) - 9  mil  Inner  gap  (10  modes) 
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Figure  4.8:  Measured  permittivity  and  permeability  of  MagRAM  are  shown  as 

er  =  e'r  —  je"  and  fir  =  fi'r  —  j/i",  repectively.  Error  bars  indicate  potential  error  in 
radial  width  (±2  mil)  and  length  (±2  mil)  measurements.  Using  the  modal  method 
with  10  modes  considered  produces  a  very  accurate  correction  for  the  real  (a)  and 
imaginary  (b)  part  of  the  permittivity  and  real  (c)  and  imaginary  (d)  part  of  the 
permeability. 
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2  mil  radial  width  measurement  error 


2  mil  sample  length  measurement  error 


(c) 


(d) 


Figure  4.9:  Measured  permittivity  and  permeability  of  MagRAM  are  shown  as 

er  =  e'r  —  je"  and  fir  =  [i!r  —  jfi",  repectively.  Error  bars  indicate  potential  error  in 
radial  width  (±2  mil)  and  length  (±2  mil)  measurements.  Using  the  modal  method 
with  10  modes  considered  produces  a  very  accurate  correction  for  the  real  (a)  and 
imaginary  (b)  part  of  the  permittivity  and  real  (c)  and  imaginary  (d)  part  of  the 
permeability. 
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Figure  4.10:  Calculated  permittivity  and  permeability  of  MagRAM  are  shown  as 
er  =  e'r  —  je"  and  nr  —  i>!r  —  j/i",  repectively.  Error  bars  indicate  potential  error  in 
sample  radial  width  (±2  mil)  and  length  (±2  mil)  measurements.  Results  are  for 
measurements  at  1  GHz  and  increasing  number  of  modes.  The  measurements  are 
reasonable,  but  appear  to  diverge  after  20  modes. 
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- Fully  Filled 


20  modes  — —  Lumped  Circuit 


Figure  4.11:  Calculated  permittivity  and  permeability  of  MagRAM  are  shown  as 
er  =  e'r  —  j e"  and  /ir  =  p!r  —  j/i",  repectively.  The  outer  air  gap  is  varied  (by  using 
multiple  samples)  from  3  mil  to  12  mil.  The  lumped  circuit  is  shown  for  comparison. 
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Applying  both  root  search  methods  to  a  complex  valued  system,  e'  —  je"  for 
example,  introduces  another  degree  of  freedom  to  the  root  search  as  both  the  real  and 
imaginary  parts  are  variable.  Both  root  search  methods  test  the  magnitude  of  the 
difference  in  sequential  steps  in  the  interpolation,  thereby  removing  the  additional 
degree  of  freedom  when  determining  convergence. 

A  complex  value  root  search  may  produce  non-physical  results  for  very  low 
loss  or  lossless  materials.  For  example,  e"  is  defined  to  always  be  positive,  thereby 
matching  the  physical  system.  The  root  search  in  this  research  was  not  constrained 
to  this  physical  expectation.  In  testing  for  convergence,  the  threshold  provided  by 
the  user  may  by  reached  when  the  imaginary  part  of  the  permittivity  is  very  small, 
but  negative.  Therefore,  very  small  (<  .05)  negative  values  of  e"  are  reported  on 
occasion.  Allowing  the  algorithm  results  to  vary  unconstrained  provides  a  good  check 
of  the  algorithm.  If  the  value  of  e",  for  example,  becomes  very  negative,  it  is  a  visual 
indication  that  either  the  model  assumptions  were  violated  or  the  algorithm  is  not 
performing  correctly. 

4-4-1  Newton  Root  Search  for  Minimization.  A  Newton  root  search  was 
used  to  perform  the  minimization  in  (4.1).  An  accurate  initial  guess  is  always  a 
concern  when  using  a  Newton  root  search.  In  this  application,  the  user  is  expected  to 
have  a  reasonable  guess  for  the  permittivity  and  permeability  of  the  material  being 
measured.  If  multiple  roots  exist  in  (4.1),  the  Newton  root  search  may  find  the  wrong 
root  if  the  initial  guess  is  not  close  enough  to  the  actual  value.  For  the  non-magnetic 
Rexolite®  material  with  a  100  mil  outer  air  gap,  the  magnitude  of  the  difference 
between  the  theoretical  and  experimental  S2i  were  evaluated  and  plotted  (figure  4.12) 
in  a  complex  region  immediately  surrounding  the  expected  permittivity  to  look  for 
potential  multiple  roots  for  both  1  mode  and  10  mode  solutions2.  No  extra  roots  were 
noted  in  either  S2i  or  Sn  plots  and  there  is  little  difference  in  the  regions  calculated 

2For  non-magnetic  material,  =  1  is  assumed.  Therefore,  only  a  one  dimension  root  search  is 
needed  using  either  S21  or  Sn- 
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using  1  mode  versus  10  modes  with  the  exception  that  the  solution  is  shifted  more 
toward  the  actual  value  (as  expected  when  increasing  the  number  of  modes).  This 
indicates  that,  for  a  non-magnetic  material,  only  a  reasonable  initial  guess  is  necessary. 

Analyzing  the  Newton  root  search’s  sensitivity  to  an  initial  guess  is  more  difficult 
when  magnetic  material  is  evaluated.  A  two  dimensional  root  search  is  used  and  hence 
the  permittivity  and  permeability  are  iterated  with  (4.1).  The  simple  plots  produced 
for  Rexolite®  are  not  possible.  Rather  an  empirical  analysis  is  more  appropriate 
and  understandable.  Initial  guesses  of  er  =  1  —  1  j  and  /zr  =  1  —  1  j  produced  the 
same  result  as  a  more  accurate  initial  guess.  It  must  be  noted  that,  for  the  sake  of 
efficiency  and  accuracy,  the  initial  guess  is  only  used  to  find  the  first  values  of  er  and 
Hr  for  the  material.  Calculations  for  the  next  higher  frequency  use  as  the  initial  guess 
the  result  from  the  previous  frequency.  This  technique  assumes  the  permittivity  and 
permeability  do  not  change  rapidly  with  frequency. 

The  Newton  root  search  is  appropriate  for  evaluating  (4.1).  The  roots  are 
distinct  and  separated.  Provided  the  initial  guess  is  reasonable,  the  Newton  algorithm 
converged  to  the  necessary  root. 

4.4.2  Muller  Root  Search  for  Wavenumbers.  Determining  the  wavenumbers 
in  Region  II  of  the  coaxial  test  fixture  is  the  primary  and  most  critical  calculation  in 
the  modal  method.  The  wavenumbers  are  determined  by  finding  the  roots  of  (3.21), 
reprinted  here  as 

kpinZin  Vo  (kPlnR)  W1  (kP2nR)  —  kP2nZ2n  1 1  ( kPlnR )  Wo  ( kP2nR )  =  0  (4-2) 

The  oscillatory  nature  of  the  Bessel  functions  guarantees  an  infinite  number  of  roots 
that  must  be  determined  in  sequential  order.  In  the  case  of  heavily  loaded  materials, 
the  first  roots  in  the  sequence  are  often  in  close  proximity  to  each  other.  For  this 
reason,  a  Muller  root  search  method  was  used  to  ensure  the  proper  root  was  located 
in  sequential  order  [18]. 
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Figure  4.12:  For  a  Rexolite®  sample  with  100  mil  outer  air  gap,  the  magnitude, 
in  dB,  of  the  difference  between  the  theoretical  and  experimentally  measured  S21 
parameter  is  plotted  as  a  complex  surface  around  the  expected  solution  of  er  =  2.525. 
The  smallest  value  on  the  surface  (blue)  is  the  solution  of  the  matrix  (and  hence  the 
permittivity  of  the  material).  A  single  root  is  apparent  for  the  modal  method  using 
1  mode  (a)  and  10  modes  (b).  Similar  results  were  noted  for  the  Sn  parameter. 
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The  Muller  root  search  method  requires  three  initial  guesses.  The  roots  of  (4.2) 
vary  based  on  material  parameters,  air  gap  size,  test  fixture  dimensions,  and  frequency. 
For  low  frequency  measurements  of  low-loss,  low  dielectric  constant  material  such  as 
Rexolite®,  the  initial  guess  of  the  roots  of  (4.2)  is  sufficiently  estimated  with  the  roots 
of  (3.20),  the  fully  filled  coaxial  line  case.  Table  2.2  validates  this  assumption.  There 
is  very  little  difference  between  the  filled  and  air  gap  wavenumbers. 

Wavenumbers  for  a  heavily  loaded  material  may  be  significantly  different  be¬ 
tween  a  filled  and  air  gap  scenario,  even  for  small  air  gaps  as  shown  in  table  2.3  and 
figure  4.13.  An  initial  guess  in  this  case  is  extremely  important.  For  heavily-loaded 
MagRAM,  the  initial  guess  was  calculated  based  on  asymptotic  expansions  for  (4.2) 
at  the  lowest  measurement  frequency,  50  MHz.  The  frequency  was  increased  incre¬ 
mentally  by  15  MHz  steps,  tracking  the  root  at  each  step,  until  the  desired  frequency 
was  reached  shown  in  figure  4.13. 

Initial  guesses  for  for  high  frequency  and  heavily  loaded  material  samples  are 
made  with  asymptotic  expansions  as  goes  to  oo.  Large  argument  approximations 
of  Bessel  functions  are  well-known  [1]  as 


Jn[z) 


7T 

2 


Vn(z) 


7 T 
2 


for  8z  »  4 n2  —  1.  Substituting  into  (4.2)  and  reducing  produces  a  large  argument 
approximation  as 


^  ^  tan  (y/j*  +  k?(R  -  a)^  + 

^  tan  (Vbn  +  k%(b  -  fyj  =  0  (4.3) 

The  tangent  function  is  not  considered  stable  based  on  its  divergence  at  multiples 
of  |,  but  that  fact  is  utilized  in  estimating  the  roots.  A  heavily  loaded  material 
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produces  a  high  contrast  between  the  air  gap  and  the  material  (i.e.  for  an  outer  air 
gap  k\  »  k2).  Additionally,  considering  the  small  air  gaps  used  in  this  research, 
the  radial  width  of  the  air  gap  is  significantly  smaller  than  the  radial  width  of  the 
material  (i.e.  for  an  outer  gap  ( b  —  R )  <<  (R  —  a)).  The  tangent  function  arguments 
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n 


(4.4) 


K  =  i/K  +  *Z(b-R) 


(4.5) 


and  the  tangent  function  weights  are  therefore  significantly  different,  are  related  via 
7„,  and  ensure  the  two  tangent  functions  are  separated  in  magnitude  and  phase.  The 
large  difference  in  magnitude  requires  that  the  roots  of  4.3  are  at  points  very  near  to 
poles  of  the  tangent  functions  such  that  large  %  is  estimated  reasonably  well  as 
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(4.7) 


where  4>n  =  |  +  mi  and  6n  =  |  +  nil.  The  gamma  estimate  must  be  based  on  the 
material  properties  and  thus  (4.6)  is  used  for  outer  air  gap  gamma  estimates  and  (4.7) 
for  inner  air  gaps. 

The  first  root  of  4.2  at  low  frequencies  for  the  heavily  loaded  MagR.AM  falls  far 
below  the  criteria  for  a  large  argument.  Therefore,  a  small  argument  approximation 
(table  F.l)  is  used  to  reduce  (4.2)  to 
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-  k\ 

1  -  Q 
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where 


n  = 


-er2  Mf) 

erl  In® 


(4.9) 


The  n  subscript  is  absent.  Only  the  principle  value  of  the  small  argument  approxi¬ 
mation  for  Yn(z)  is  used  (where  z  is  complex).  The  small  argument  approximation 
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provides  an  initial  guess  for  the  first  zero  of  (4.2)  at  low  frequency  with  a  heavily 
loaded  material  and  therefore  multiple  values  are  not  needed. 

The  initial  guesses  for  %  calculated  from  asymptotic  techniques,  in  conjunction 
with  the  Muller  root  search  method,  are  able  to  quickly  and  accurately  find  the 
actual  value  of  jn.  Figure  4.14  is  a  complex  surface  of  values  calculated  via  (4.2)  for 
varying  values  of  %  at  2GHz  for  a  coaxial  line  fully  filled  with  MagRAM  (a)  and 
a  coaxial  line  with  MagRAM  having  a  9  mil  inner  air  gap  (b).  The  inner  air  gap 
case  is  used  because  it  proved  to  be  the  most  difficult  to  track  over  frequency.  The 
first  four  wavenumbers  determined  by  the  Muller  method  and  the  complex  path  the 
Muller  method  took  to  find  the  root  are  both  shown.  The  lowest  magnitudes  in  the 
plot  (blues  spots)  are  the  actual  roots  the  Muller  method  is  attempting  to  locate. 
Although  the  roots  are  significantly  different  from  the  no  gap  to  9  mil  inner  air  gap 
case,  with  the  proper  initial  guess  the  Muller  method  quickly  locates  the  roots.  This 
is  true  despite  the  decrease  in  gradient  between  successive  roots  in  the  9  mil  inner  air 
gap  case  in  comparison  to  the  no  air  gap  case. 

The  accuracy  of  the  Bessel  function  calculations  are  critical  in  determining  77 
from  (4.2).  The  Matlab®  software  package  was  used  to  implement  all  of  the  routines 
required  for  this  research.  The  native  Bessel  functions  were  used.  The  first  50  wave 
numbers  for  a  MagRAM  material  with  9  mil  inner  air  gap  at  50MHz  were  calculated 
and  are  plotted  in  figure  4.15.  The  initial  guess,  7®st,  is  also  plotted,  ft  appears 
that  (4.2)  is  being  approximated  by  a  tangent  function  for  the  higher  order  modes 
(large  argument).  This  is  consistent  with  the  large  argument  approximation  used  in 
the  derivation  of  initial  guesses,  implying  that  the  large  argument  approximations  are 
being  used  in  the  software  package  when  calculating  Bessel  functions.  Based  on  the 
figure,  including  more  than  23  modes  is  counter  productive  and,  not  knowing  at  what 
point  the  argument  is  large  enough  that  the  software  package  begins  to  approximate 
the  Bessel  function,  including  more  than  15  modes  may  alter  the  outcome  for  complex 
7 n.  This  is  evident  in  the  modal  method  results  for  MagRAM  (figure  4.10),  but  not 
in  the  results  for  Rexolite®.  This  implies  that  the  approximation  is  embedded  in 
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Figure  4.13:  The  real  (a)  and  imaginary  (b)  parts  of  the  wavenumber,  7n,  for  the 
first  four  modes,  TM0n,  as  a  function  of  frequency.  The  initial  guesses,  estimated 
using  asymptotic  techniques,  are  shown.  For  70,  the  initial  guess  begins  to  diverge 
from  the  actual  value  at  500  MHz. 
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Figure  4.14:  The  potential  solutions  to  the  characteristic  equations  for  a  fully  filled 
(a)  and  air  gap  (b)  coaxial  line  plotted  as  a  complex  plane.  The  Muller  root  search 
locates  the  lowest  point  on  the  surface  (the  blue  areas).  The  frequency  is  2  GHz 
for  these  plots.  When  the  MagRAM  fully  fills  the  coaxial  line  the  wavenumbers  are 
visually  sequential  and,  with  an  accurate  initial  guess,  are  found  quickly.  Conversely, 
when  an  air  gap  is  present,  is  not  necessarily  sequential  and  the  gradient  between 
successive  wavenumbers  is  shallower,  implying  a  more  accurate  initial  guess  is  critical. 
The  first  four  wavenumbers  and  the  search  path  followed  by  the  Muller  root  finding 
method  are  labeled  on  the  graphs. 
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the  native  Bessel  functions  when  the  large  argument  is  complex  (the  MagRAM  case), 
but  not  for  purely  real  (as  in  the  Rexolite®  case).  This  requires  further  research  to 
validate. 
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Figure  4.15:  The  real  (a)  and  imaginary  (b)  parts  of  the  wavenumbers  ,  7„,  for 

n  =  0, 2, . . . ,  49  at  50MHz.  Wavenumbers  >  23  are  non-physical  in  that  the  imaginary 
part  is  <  0.  Note  the  initial  guesses,  7®st,  remain  physical  at  higher  order  modes,  but 
the  located  roots  appear  to  follow  a  tangent  function. 
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V.  Conclusions  and  Recommendations  for  Improvement 

The  excitation  of  higher  order  modes  is  a  major  source  of  error  in  material  charac¬ 
terization  measurements  using  the  NRW  algorithm.  Higher  order  modes  in  coaxial 
line  measurements  are  typically  caused  by  axially  symmetric  air  gaps  between  the 
material  sample  and  the  inner  or  outer  conductor.  A  modal  method  was  developed 
assuming  axially  symmetric  modes  are  excited.  The  method  was  proven  to  remove  er¬ 
ror  associated  only  with  axially  symmetric  modes.  One  property  of  the  modal  method 
is  the  convergence  of  the  result  to  the  actual  material  properteris  as  the  number  of 
modes  are  increased,  provided  the  only  significant  measurement  error  is  from  axially 
symmetric  higher  order  mode  excitation.  This  is  in  contrast  to  the  lumped  circuit 
correction  factor  that  weights  the  NRW  result  based  on  material  and  air  gap  dimen¬ 
sions.  The  weighting,  although  fairly  accurate,  does  not  provide  confidence  in  the 
solution,  particularly  for  heavily  loaded  materials. 

5.1  Recommended  Improvements 

The  shortcomings  of  the  modal  method  developed  here  and  potential  ways  to 
account  for  them  are  described  here. 

5.1.1  Initial  Guesses  for  Wavenumbers.  Using  large  and  small  argument 
approximations  for  Bessel  functions  proved  very  valuable  for  estimating  the  wavenum¬ 
bers  in  heavily  loaded  MagRAM  with  air  gaps  considered  in  this  research.  The  method 
is  simplified  at  best  and  would  benefit  greatly  from  a  more  thorough  exploration  and 
mathematical  analysis.  The  wavenumbers  are  extremely  important  and  as  such  the 
initial  guesses  are  a  critical  component  of  the  modal  method. 

5.1.2  Compensation  for  Eccentric  Center  Conductor.  A  properly  machined 
sample  that  fully  fills  the  cross-section  of  the  coaxial  line  (i.e.  no  air  gap)  physically 
supports  the  center  conductor  in  the  center  of  the  waveguide.  As  shown,  if  the  sample, 
when  measured,  has  either  an  inner  or  outer  air  gap  the  physical  support  that  keeps 
the  center  conductor  on  the  center  axis  of  the  coaxial  line  is  gone  (see  figure  5.1).  In 
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Figure  5.1:  The  physical  support  provided  by  a  sample  that  fully  fills  the  coaxial 
line  is  lost  when  an  air  gap  is  present.  More  physical  support  to  the  conductor,  a 
different  mathematical  model,  or  a  combination  are  needed  to  overcome  the  effect  on 
measurements. 

the  case  of  an  outer  air  gap,  not  only  is  the  outer  support  removed,  but  the  sample  puts 
additional  weight  on  the  center  conductor.  This  research  assumed  a  center  conductor 
that  remained  along  the  center  axis  when  an  air  gap  was  present  and  therefore  only 
axially  symmetric  modes  were  excited.  Additional  supports  are  necessary  within  the 
coaxial  line  to  ensure  the  center  conductor  remained  in  the  center.  However,  given  the 
dimensions  in  figure  5.2,  for  a  small  droop,  d,  a  perturbation  analysis  of  the  fields  in 
the  material  and  air  gap  region  could  be  performed  to  produce  a  modified  propagation 
constant,  7n,  for  use  in  the  modal  algorithm  defined  in  this  research. 

5.1.3  Compensation  for  Eccentric  Material.  The  presence  of  an  inner  and 
outer  air  gap  in  the  material  region  was  not  addressed  by  this  research.  Assuming 
the  coaxial  center  conductor  does  not  droop,  an  inner  and  outer  air  gap  will  produce 
the  scenario  shown  in  figure  5.3.  Note  that  for  a  small  inner  radius  of  the  material, 
ri,  the  droop,  d,  will  be  small  and  a  perturbation  analysis  similar  to  the  previous 
section  could  be  accomplished.  However,  as  d  increases  the  held  modes  are  less  axially 
symmetric.  For  a  significant  d ,  a  rigorous  means  of  determining  is  found  via 
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Figure  5.2:  An  air  gap  in  the  measurement  reduces  the  physical  support  to  the 

center  conductor.  Higher  order  modes  that  are  not  axially  symmetric  may  be  excited. 
The  dimensions,  particularly  the  droop,  d,  are  exaggerated  in  the  drawing  for  clarity. 

the  addition  theorem  of  Bessel  functions.  The  addition  theorem  translates  the  field 
definitions  from  the  origin  of  the  material  to  the  origin  of  the  conductors  (or  vice 
versa).  The  application  of  the  addition  theorem  involves  another  truncated  series  and 
significantly  more  terms  to  calculate.  This  method  has  previously  been  applied  to 
determining  the  wavenumbers  for  eccentric  center  conductors  [13],  [19]  and  eccentric 
circular  dielectric  rods  [20].  With  this  method,  a  single  sample  could  be  machined  that 
could  be  measured  at  room  temperature  as  well  as  higher  temperatures  (provided  the 
conductors  did  not  expand  to  the  point  of  applying  excessive  pressure  on  the  sample). 
Beyond  the  mathematical  difficulty  in  applying  the  addition  theorem,  there  is  also 
the  physical  difficulty  of  measuring  the  dimensions  ry,  r2,  a,  b ,  and  d  as  the  sample 
and  conductor  are  heated.  A  method  must  be  devised  to  accurately  measure  these 
dimensions  during  the  testing  process  or  a  good  estimate  provided. 

5.1.4  Impedance  Boundary  Conditions.  The  coaxial  test  fixture  inner  and 
outer  conductors  are  assumed  to  be  PEC  (er  — >  00).  Considering  the  electric  field  will 
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Figure  5.3:  An  air  gap  in  the  measurement  reduces  the  physical  support  to  the 

center  conductor.  Higher  order  modes  that  are  not  axially  symmetric  may  be  excited. 
The  dimensions,  particularly  the  droop,  d,  are  exaggerated  in  the  drawing  for  clarity. 

penetrate  the  conducting  walls  like 


2 


cu/xcr 


when  a  is  finite,  but  large,  and  the  frequency  ,u>,  is  large,  the  electric  field  within 
the  conductors  are  easily  approximated  as  zero.  However,  as  the  frequency  decreases 
and  a  remains  finite,  the  electric  field  will  begin  to  penetrate  the  conducting  walls 
and  power  will  be  dissipated  in  the  walls  based  on  S  =  E  x  H*.  The  tangential 
electric  field,  E2  for  a  TMgn  mode,  must  remain  continuous  across  the  boundary, 
h  x  Ei  =  h  x  E2,  and  therefore  the  electric  field  penetration  of  the  conducting  wall 
and  the  subsequent  attenuation  has  a  direct  impact  on  the  intensity  of  the  fields 
within  the  waveguide.  The  effect  is  seen  by  making  measurements  without  a  sample 
present.  The  result  is  expected  to  be  equivalent  to  the  properties  of  freespace,  but 
figure  5.4  shows  the  loss  of  the  coaxial  line  is  more  noticeable  at  lower  frequencies.  This 
impacts  the  accuracy  of  the  measurements  as  evident  in  the  low  frequency  regimes 
of  the  data  presented  in  this  research.  The  effect  is  expected  to  be  more  pronounced 
as  the  loading  of  the  material  increases.  Reformulating  the  results  of  this  research 
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Figure  5.4:  The  measured  real  (a)  and  imaginary  (b)  part  of  the  permittivity  and 
permeability  for  the  coaxial  test  line  when  no  sample  is  present.  The  result  should  be 
that  of  freespace  (er  =  1,  pr  =  1).  The  low  frequency  perturbations  are  likely  due  to 
the  finite  impedance  of  the  coaxial  line. 

using  impedance  boundary  conditions  is  expected  to  improve  the  accuracy  of  lower 
frequency  measurements  in  the  coaxial  test  fixture. 

5.1.5  Time  Domain  Analysis  of  Scattering  Parameters.  The  generally  large 
difference  between  sequential  wavenumbers  provides  an  excellent  opportunity  to  an¬ 
alyze  the  scattering  of  individual  modes.  In  time,  each  mode  will  propagate  at  a 
different  rate  based  on  the  group  velocity,  as  a  function  of  the  wavenumber.  Sampling 
the  fields  in  time  or  increasing  the  frequency  bandwidth  of  the  measurement  should  al¬ 
low  the  different  propagating  modes  to  be  identifiable  such  that  they  could  be  treated 
individually  in  the  modal  matrix.  The  first  non-axially  symmetric  TM  mode  is  TM^q 
with  a  cutoff  frequency  at  14  GHz  in  a  free  space  filled  coax  line  (see  table  2.1).  This 
offers  a  potential  four  fold  increase  in  measurement  bandwidth  (assuming  a  pure  TEM 
excitation). 

5.1.6  Non-linear  Least  Squares  Algorithm  for  Minimization.  The  Newton 
root  search  is  performed  well  for  this  research,  but  is  not  the  only  algorithm  available 
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to  accomplish  the  minimization  of  (3.1).  Previous  research  has  used  a  non-linear 
least  squares  algorithm  to  minimize  the  difference  between  theoretical  and  measured 
scattering  parameters.  One  advantage  is  the  use  of  all  four  scattering  parameters, 
S“p,  S^15,  S^,  and  S^,  in  the  process.  For  a  simple  material  and  good  forward  and 
reverse  measurements,  the  inclusion  of  all  four  scattering  parameters  is  expected  to 
stabilize  the  results. 

5.1.7  Efficiency  in  calculation  of  Bessel  Functions.  Matlab®  is  not  de¬ 
signed  for  computational  efficiency  and  the  computation  of  Bessel  functions  are  no 
exception.  Two  dimension  root  searches  using  10  modes,  particularly  for  magnetic 
materials,  easily  took  5-6  hours  to  evaluate.  The  primary  calculations  responsible  for 
the  amount  of  time  are  the  Bessel  functions.  For  this  research  to  be  practical  in  a 
testing  environment,  the  code  must  be  ported  to  a  language  that  supports  very  fast 
calculation  of  Bessel  functions  (i.e.  Fortran®). 
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Appendix  A.  Wave  Equation 

Maxwell’s  Equations  are  a  collection  of  experimentally  derived,  mathematical  descrip¬ 
tions  for  electromagnetic  waves.  In  time-harmonic  (frequency  domain)  form  they  are 
given  as 


V  x  E  =  -jcuB 

Faraday’s  Law 

(A.l) 

V  x  H  =  J  +  jo;D 

Ampere’s  Law 

(A.2) 

V  •  D  =  qe 

Gauss’s  Law  for  D 

(A.3) 

V  ■  B  =  0 

Gauss’s  Law  for  B 

(A.4) 

where  qe  is  the  static  charge  in  the  system  and  u ;  is  the  radian  frequency  of  the  held. 

In  words,  Faraday’s  law  states  that  an  electric  held,  E.  rotates  around  a  mag¬ 
netic  hux  density,  B.  Ampere’s  law  states  that  a  magnetic  held,  H.  rotates  around  a 
linear  combination  of  an  electric  current,  J,  and  an  electric  hux  density,  D.  Gauss’s 
Law  states  that  the  electric  hux  density,  D,  diverges  from  an  electric  charge  and  the 
magnetic  hux  density,  B,  does  not  diverge  from  a  magnetic  charge,  implying  that 
the  magnetic  held  is  purely  rotational,  divergence  of  a  held  integrated  over  a  closed 
surface  is  equal  to  the  electric  or  magnetic  charge  contained  within  that  surface  and 
is  applicable  to  both  the  electric  and  magnetic  hux  densities  respectively. 

The  test  setup  of  network  analyzer  and  coaxial  line  test  hxture  are  assumed  free 
of  static  sources.  The  excitation  is  generated  external  to  the  coaxial  line  and  therefore 

oii  J*  =  0. 

Further,  it  is  assumed  that  the  material  is  linear,  homogeneous,  and  isotropic 
and  therefore  the  electric  and  magnetic  hux  densities  are  a  function  of  the  material 
in  which  the  respective  helds  are  applied.  The  following  constitutive  relationships 


D  =  eE 

(A.5) 

B  =  pH 

(A.6) 
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are  therefore  valid  within  the  material  where  e  is  the  permittivity  and  /i  the  perme¬ 
ability  of  the  material. 

Maxwell’s  equations  can  now  be  simplified  using  this  assumption.  Removing 
the  sources  and  substituting  (A. 5)  into  (A.l)  and  (A. 3),  and  (A. 6)  into  (A. 2)  and 
(A. 4)  yields 


V  x  E  =  —juj/jH. 

(A.7) 

V  x  H  =  jcoe  E 

(A.8) 

V  •  E  =  0 

(A.9) 

V  •  H  =  0 

(A. 10) 

or  specifically  a  system  of  coupled,  first  order  partial  differential  equations. 

Although  (A. 7)- (A.  10)  completely  describe  the  fields  with  the  coaxial  line,  the 
mathematical  solution  is  difficult  to  determine.  The  system  of  equations  are  simplified 
and  uncoupled  by  application  of  vector  identities.  Applying  the  curl  operation  to 
Faraday’s  law  produces 

V  x  (V  x  E)  =  V  x  (—ju/jH.) 

=  'V(-jun)  x  H  +  (— jun)S7  x  H  (A. 11) 

=  — jcj/u(VxH)  (A. 12) 

where  the  RHS  is  expanded  via  vector  identity  [3] 

V  x  (V’A)  =  x  A  +  t/’V  x  A 

and  reduce  by  recognizing  that  the  gradient  of  a  constant  is  equal  to  zero.  This  form 
allows  substituting  (A. 8)  into  (A.  12).  Additionally  using  vector  identity 


V  x  (V  x  A)  =  V(V  ■  A)  -  V2A 


to  expand  the  LHS  of  (A.  12)  and  substituting  (A. 9)  produces 


V(0)  —  V2E  =  —joJfJi  (jcoeEi)  . 


(A.13) 


Substituting  k  =  u;2e/x  and  rearranging  into  the  standard  form  of  the  wave  equation 


V2E  +  k2  E  =  0 


(A.14) 


is  an  uncoupled,  homogeneous,  second-order  partial  differential  equation  that  de¬ 
scribes  the  behavior  of  an  electric  held  in  general.  A  similar  procedure  is  accomplished 
to  derive  the  wave  equation  as  a  function  of  the  magnetic  held 


V2H  +  fc2H  =  0 


(A. 15) 
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Appendix  B.  Solution  to  Helmholtz’s  Wave  Equation  in 
Cylindrical  Coordinates 

The  coaxial  line  is  naturally  described  in  cylindrical  coordinates.  A  solution  to  Hcl- 
moholtz’s  wave  equation  [9] 

V2T  +  &2T  =  0  (B.l) 


for  cylindrical  coordinates  is  needed  to  describe  the  fields  in  the  coaxial  line.  Expanded 
in  cylindrical  coordinates,  (B.l)  is 


1  d  (  <9T\ 
P  dp  V  dp  ) 


1  <92T  <92T  r9 

+  7W  +  ^  + 


(B.2) 


Assume  that  the  solution  T  is  separable  into  a  product  of  functions  that  are  dependent 
on  a  single  coordinate  direction  (the  well-known  separation  of  variables  technique  [1]) 
and  can  therefore  be  written  as 


*{p,<l>,z)  =  F{p)G(<l>)H{z)  (B.3) 


Using  (B.3)  in  (B.2)  produces 


1  d  (  dF\  1  / 1  d2G\  {  l  d2H\  9  „ 

¥pif>  (flp )  +  p  )  +  \H~d? )  + 


(B.4) 


The  second  and  third  parenthetical  terms  are  dependent  only  on  (j)  and  z  respectively. 
For  the  solution  to  sum  to  zero  over  all  space,  these  terms  must  also  be  constant  [9]. 
Therefore,  the  assignments 


1  d2G 
~G~dxjE 

1  d2H 
H~dW 


-m 


(B.5) 

(B.6) 


are  made  and  substituted  into  (B.4)  resulting  in 


(B.7) 


Defining 


k2p  =  k 2  -  k\  (B.8) 

and  rearranging  terms,  the  defining  equations  (B.5),  (B.6),  and  (B.7)  for  G ,  H,  and 
F  respectively  are 


%  ) +  ~  "d  F 


d2G 

d(j): 


dz2 


—  m 2]  F 

=  0 

(B.9) 

j  +  Gm2 

=  0 

(B.10) 

1 

2  +  Hkl 

=  0 

(B.ll) 

of  which  the  solutions  are  constrained  by  (B.8). 

Equation  (B.9)  is  recognizable  as  the  Bessel  differential  equation.  The  general 
solution  will  be  a  linear  combination  of  Bessel  functions  chosen  to  best  describe  the 
p  directed  held  in  the  original  problem.  For  a  coaxial  line,  the  held  will  be  a  standing 
wave  and  therefore  the  general  solution  is  defined  as 


F(p)  =  AJm  ( kpp )  +  BYm  ( kpp )  (B.12) 

where  Jm  is  a  Bessel  function  of  the  hrst  kind  with  order  m  and  Ym  is  a  Bessel  function 
of  the  second  kind  with  order  m  [9]. 

Solutions  of  the  ordinary  differential  equations  (B.10)  and  (B.ll)  are  well- 
known.  The  general  solutions  are  chosen  based  on  the  description  of  the  original 
problem.  In  the  case  of  a  coaxial  waveguide  whose  center  axis  lies  along  the  z  axis, 
the  natural  solutions  are 


G(d>) 

H(z) 


C  cos (m(j))  +  D  sin (mcj)) 
Ee~jk*z  +  FejkzZ 
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(B.13) 

(B.14) 


where  z  dependent  exponential  terms  describe  propagation  in  the  z  direction  and  (j) 
dependent  terms  describe  the  circumferential  variation  of  the  held  around  the  z  axis. 

The  solution  to  (B.l)  is  now  given  by  a  product  of  the  individual  solutions  as 
in  (B.3), 

^ m,n  (P)  0;  *0  [An,n'4n  (^pP)  T  (^pP)]  [Cm,ra  COS  T  Dm  n  sill  X 

[Em,ne-lz  +  Fm,ne*z]  (B.15) 

where  7  =  jkz  and  subscripts  m  and  n  are  added  for  denoting  unique  solutions. 
Additionally  the  constraint  equation 

k2p  =  k2  +  q2  (B.16) 


must  hold  for  all  solutions. 
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Appendix  C.  Cutoff  Frequency  of  a  TEM  Field 

The  cutoff  frequency  of  a  TEM  field  is  0  Hz.  This  property  allows  low  frequency 
measurements  when  a  TEM  field  is  excited  (such  as  in  a  two  conductor  transmission 
line).  Proving  the  cutoff  frequency  bound  is  done  by  examining  the  fields  with  respect 
to  Maxwell’s  equations,  and,  by  implication,  the  wave  equation.  By  definition,  the  z 
directed  field  components  of  a  TEM21  wave  are  equivalent  to  zero.  The  electric  and 
magnetic  fields  in  a  TEM2  wave  are  therefore  completely  described  by  transverse  field 
components 


E  =  E  t(p,<f>,z)  (C.l) 

H  =  H  t(p,(f>,z)  (C.2) 

which  must  satisfy  Maxwell’s  equations  and,  as  a  consequence,  the  wave  equation 
(A.  14). 

Faraday’s  law  is  reduced  to 


V  x  Et  =  0  (C.3) 

when  only  the  transverse  field  components  are  present.  Using  the  vector  identitiy  [3]: 

V  x  (VtVO  =  o  (C.4) 

the  transverse  electric  field  can  be  posed  as  the  transverse  gradient  of  the  scalar 
ijj(p,(j),z).  Substituting  into  Gauss’  law  and  assuming  a  simple  media 

V-Et  =  V-(Vt^) 

=  Vfy 

=  0  (C.5) 
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where  V*  is  the  gradient  of  the  transverse  components  only  or 


V 


„  d  »1  9  d 

^  dp  pdcj)^  dz 
vt 


The  wave  equation  where  only  transverse  fields  are  present  is 


V2Et  +  k2Et  =  0 


(C.6) 


(C.7) 


For  a  wave  travelling  in  the  I  direction,  the  tangential  electric  held  Et(p,(/),z)  is 
equivalent  to  Et(p,  cj))e~lz.  The  vector  Laplacian  defined  with  respect  to  transverse 
and  longitudinal  components  is 


V2  =  V2  + 


cf_ 

dz2 


(C.8) 


and  when  applied  to  the  electric  held  is 

V2E =  fv2+|^)  Etfafle-v 

=  V2Et(p,0)e-^  +  7E  t{pA)e-lz 

=  (V2  +  72)  Et(p,  0,  z)  (C.9) 


Substituting  (C.9)  into  (C.7)  results  in 

V2Ef  +  (y2  +  k2)Et  =  0 


(C.10) 


Restating  previous  results  from  (C.5) 


Ef  =  \7ti> 
=  0 
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and  substituting  into  (C.10)  produces 


V2  (VWO  +  (72  +  k2)  (V^)  =  0 
Vt  (VtV  +  (72  +  =  0 

(7  2  +  k2)Wt^  =  0 
(72  +  fc2)Ei  =  0 

The  non-trivial  result  requires  that  7s  =  —k2.  Recalling  the  constraint  equation 
(B.16)  and  k  =  it  is  evident  kp  =  0  for  a  TEM  wave.  As  such  the  propagation 

constant,  7,  varies  directly  with  frequency,  u,  and  there  is  no  frequency  at  which  the 
term  describing  propagation,  e-72,  is  purely  real.  A  purely  real  result  would  produces 
an  evanescent  (decaying)  wave.  A  TEM  wave,  when  possible,  will  always  propagate 
without  respect  to  frequency. 
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Appendix  D.  TEM  and  TM  Mode  Orthogonality 

Mode  orthogonality  is  a  critical  component  to  understanding  the  power  carrying  mech¬ 
anisms  in  waveguide.  Each  mode  must  independently  carry  power  to  which  the  total 
power  is  the  summation  of  the  individual  mode  powers  [6].  Therefore,  each  mode  func¬ 
tion  must  be  orthogonal.  Mathematically,  the  solutions  to  the  scalar  wave  equation 
must  satisfy  [1] 

l^gyhdS  =  0  ...  g  A  h  (D.l) 

Js 

or  the  equivalent  for  vectors 


VtTg  ■  Wt^hdS  =  0 


9  A  h 


(D.2) 


The  following  analysis  assumes  a  coaxial  waveguide  with  no  discontinuities  and 
only  TEM  and  TM  modes  present.  Although  a  parallel  analysis  for  TE  modes  is 
possible,  it  is  unnecessary  for  this  research. 


D.l  General  Field  Solution  in  Cylindrical  Coordinates 

Constructing  field  equations  from  vector  potentials  is  thoroughly  covered  in  [3] . 
Generating  TM2  fields  requires  only  a  z  directed  potential  function  that  satisfies  the 
scalar  wave  equation 

V2Aj  +  k2Az  =  0 

This  is  equivalent  to  (B.l)  and  therefore  (B.15)  is  a  valid  form  of  Az  or 

Az  =  %{p,(f>,z)  =  %(Pi<f>)e^z  (D.3) 

where  m,  n  pairs  are  replaced  by  a  single  mode  identifier,  g,  and  only  the  forward 
propagating  wave  is  considered  for  simplicity.  The  electric  held  for  a  TM2  mode  is 
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therefore  defined  as 


E  =  — 1 —  (k2  +  yy-)z<HJp,^)e-<z 

jweju 

=  —  [J(*:2+72)*!,(P,^)e_7'-7e-1'V1>t(,(p,^)]  (D.4) 

juicy 

where  V  is  replaced  with  the  transverse  and  longitudinal  operator  (C.6).  The  result 
is  the  separation  of  the  z  directed  field  from  the  transverse  fields.  Both  must  satisfy 
the  scalar  wave  equation. 

D.2  Orthogonality  of  z  directed  fields 

The  z  directed  fields  are  a  function  of  and  proportional  to  the  original  solution, 
(p,  e~lz .  Considering  two  unique  solutions,  (p,  cj),  z)  and  dR  (p,  0,  z),  they 
both  must  independently  satisfy  the  scalar  wave  equation  as 

V?<I>S  +  (A:2+72)>I>«,  =  0  (D.5) 

V?^+(i2  +  7l)*i  =  0 

Multiplying  each  scalar  wave  equation  by  the  opposing  solution  and  subtracting  the 
two  results  produces 


1 

juey 


(732  -7S>iyn  =  -  s'ftV?*,  (d.6) 

Thought  the  LHS  is  well-suited  to  analysis  via  (D.l),  the  RHS  can  be  further  manip¬ 
ulated  via  the  identities 


v2^  =  v-w> 


(D.7) 


(V  •  A)  =  V  •  (V’A)  -  A  • 


(D.8) 
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Coax  Cross  Section 


f/{}dS 


f  0dl 


s 


c 


Figure  D.l:  In  two  dimensions,  the  surface  and  closed  line  integration  for  the 

orthogonality  condition  in  an  empty  coaxial  is  the  freespace  region  bounded  by  the 
inner  and  outer  conductor.  The  closed  line  integration  connects  the  inner  and  outer 
conductor  via  a  twice  followed  path  whose  integration  is  equal  and  opposite  thereby 
reducing  the  the  connecting  integration  to  zero. 

to  produce 


TgVf2Th  -  VhV2tVg  =  V*  •  -  Vt  •  {Vh'VtVg) 

=  Vt  •  {*gVt*h  -  VhVtVg) 


Substituting  this  result  into  (D.6),  integrating  with  respect  to  the  paths  in  figure  D.l, 
and  applying  the  divergence  theorem  in  two  dimensions  [9] 


(792  lh)  J  J 


Vg^hdS 


s 


s 


(D.9) 


c 


Recalling  that  and  are  proportional  to  the  z  directed  electric  field  for  mode 
g  and  h  respectively,  the  enforcement  of  tangential  boundary  conditions  on  the  PEC 
wall  of  the  coaxial  line  requires  Ez  =  0  and  therefore  requiring  the  solutions  and 
xI>h  are  equivalent  to  0  along  the  integration  contour.  For  this  reason  (D.9)  reduces 
to  (D.l)  and  the  z  directed  fields  are  proven  orthogonal. 

It  is  important  to  note  the  following: 


1.  The  development  above  applies  to  TM  fields.  By  definition,  Ez  =  0  for  a  TEM 
field  and  (D.l)  reduces  to  a  trivial  result  ensuring  orthogonality. 

2.  Evident  in  (D.9)  is  the  potential  for  a  trivial  solution  should  two  unique  solu¬ 
tions  have  identical  eigenvalues  (7^  =  7^).  The  solutions  are  then  considered 
degenerate  [6].  In  assuming  only  axially  symmetric  modes,  this  research  does 
not  require  consideration  of  degenerate  modes. 

3.  The  solutions  and  dT  are  only  valid  over  the  surface  shown  in  figure  D.l. 

D.3  Orthogonality  of  transverse  fields 

The  transverse  fields  are  a  function  of  and  proportional  to  the  transverse  gra¬ 
dient  of  the  original  solution,  (p,  </>)  e-72.  Unique  solutions  are  then  (p,  (j),  z ) 
and  VjdU  (p,  </>,  z).  A  simple  derivation  is  produced  by  substituting  the  transverse 
solutions  into  (D.7)  and  (D.8)  as 


VtT9  ■  VfdU  =  Vt  ■  (dUVfTg)  -  ■  VfTg) 


=  Vt  •  (T,VfTg)  -  (D.10) 


The  integral  of  (D.10)  over  the  cross-sectional  surface  of  the  waveguide  produces 
Green’s  identity  in  two  dimensions  [9] 


V,T0  ■  Vt^hdS  =  &  •  hdl 


t^gdS 


C 


(D.ll) 


Again,  the  RHS  closed  contour  integration  is  equivalent  to  zero  based  on  the  enforce¬ 
ment  of  tangential  boundary  conditions  and  the  proportionality  of  to  Ez.  By 
substituting  (D.5)  into  (D.ll) 


IJ  VtTg  ■  Vt*hdS  =  -( k 2  +  752)  I j  VhVgdS  (D.12) 

s  s 
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for  which  the  already  proven  (D.l)  reduces  the  RHS  to  0,  proving  (D.2)  and  conse¬ 
quently  the  orthogonality  of  transverse  fields  of  a  TM2  wave.  Recalling  the  use  of 
E tem  =  Vf'h  in  appendix  C,  it  is  apparent  that  the  orthogonality  proven  for  the 
transverse  fields  of  a  TM*  wave  applies  to  a  TEM~  wave. 


Appendix  E.  Bessel  Integration  Identities 

The  following  identities  are  used  to  produce  reduced  analytical  expressions  for  the 
coupling  integrals  of  the  MAM.  The  identities  are  developed  in  similar  forms  in  [8], 
[11],  [14],  and  [16], 


/'  x2 

I  [Zn(ax)]2xdx  =  y  {[Zn(ax)]2  -  Zn_i(ax)Zn+i(ax)}  (E.l) 

x2  f  2n  1 

=  —  <  [Zn(ax)]2  +  [Zn_i(ax)]2 - Zn„i(ax)Zn(ax)  \ 

2  y  ax  J 

where  Zn  is  a  Bessel  function  or  linear  summation  of  Bessel  functions  of  the  first  and 
second  kind. 


Proof 

Starting  with  Bessel’s  differential  equation 


x2u"  +  xu'  +  [(ax)2  —  n2]  u 

=  0 

2 x2u'u"  +  2 xu'u'  +  2 uu'  (ax)2  —  2 uu'n2 

=  0  . 

. .  (mult,  by  2 v!) 

2 x2u'u"  +  2 xu'u'  +  2 uu'  (ax)2  —  2uu'n 2  +  (2 a2xu2  —  2 a2xu2) 

=  0  . 

. .  (adding  zero) 

{x2[u')2  +  [(ax)2  —  n2]  u2}  —  2 a2xu2 
ax 

=  0 

|x2(«')2  +  [(ax)2  —  n2]  u 2j 

=  2“2J 

[  xu2dx 

The  solution  to  Bessel’s  equation,  u,  must  take  the  form  of  a  Bessel  function  or 
linear  summation  of  Bessel  functions.  This  research  uses  linear  summation’s  of  Bessel 
functions  of  the  first  and  second  kind  (  Jn  and  Yn  ,  respectively)  such  that 

u  =  Zn  ( ax ) 

where  Zn  can  take  the  form  Jn  (Ax),  Yn  (Ax),  or  AJn  (Ax)  +  BYn  (Ax)  It  is  important 
to  note  that  derivatives  in  the  Bessel  differential  equation  are  taken  with  respect  to 
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x  such  as  «'  =  jh  Applying  this  to  the  substitution  u  =  Zn  (ax)  produces 


du 

dx 


—  Zn  (ax) 
dx 

a—Zn  (ax) 
dx 

aZ'n  (ax) 


Continuing  the  proof  by  substituting  for  u 


x 


[Zn  (ax)]2  xdx  =  —  ^  [Z'n  (ctx)]  +  [Zn  (ax)Y 


(^-\2[Zn(ax)f 

\ax/ 


the  recursion  relations  [1] 


Zn_x  (ax)  =  —Zn  (ax)  +  d  Zn  (ax) 

ax  d(ax) 

Zn+ 1  (ax)  =  —Zn  (ax)  -  d  Zn  (ax) 

ax  d(ax) 


are  readily  used  to  complete  the  proof  as 


[Zn (ax)} 2  xdx  =  2—  { [Zn(ax)]2  -  Zn_1(ax)Zn+1(ax)j 
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Zn(a:x)Wn(/3x)xdx 


(E.2) 


x/?Zn(ax)Wn_i(/3x)  —  xaZn_i(o;x)Wn(/3x) 

a2  —  /32 


where  Zn  and  Wn  are  Bessel  functions  or  linear  summation  of  Bessel  functions  of  the 
first  and  second  kind. 


Proof 


Starting  with  Bessel’s  differential  equation 

x2u"  +  xu'  +  [(ax)2  —  n2]  u  =  0 
x2v"  +  xv'  +  [(/ 3x )2  —  n2]  v  =  0 

where  u  and  v  are  unique  solutions.  Multiplying  by  ^  and  ^  respectively  and  sub¬ 
tracting  the  result  gives 

xu"v  +  u'v  +  (a2x2  —  n2)—  —  xv"u  —  v'u  —  ((32x2  —  n2)—  =  0 


x 


x 


xu"v  —  xv"u  +  u'v  —  uv'  +  (a2  —  /32)uvx  =  0 

{x(u'v  —  uv')\  +  uvx(a2  —  (32)  =  0 
ax 


Integrating  with  respect  to  x  produces 


The  solutions  to  Bessel’s  equation,  u  and  v,  must  take  the  form  of  a  Bessel  function  or 
linear  summation  of  Bessel  functions.  This  research  uses  linear  summation’s  of  Bessel 
functions  of  the  first  and  second  kind  (  Jn  and  Yn  ,  respectively)  such  that 


u 


v 


Zn  (ax) 
Wn  (fix) 
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where  both  Zn  and  Wn  can  take  the  form  Jn  (Ax),  Yn  (Ax),  or  AJn  (Ax)  +  BYn  (Ax). 
Using  the  recursion  relation  [1] 


—  Zn(ax )  =  aZ'  (ax)  =  a 
ax 


n 


Zn- 1  (ax) - Zn  (ax) 

ax 


and  substituting  u  =  Zn(ax)  and  v  =  Wn(/3x)  produces 


Zn(ax)Wn(f3x)xdx  = 


x  {Zn(ax)£Wn([3x)  -  £Zn(ax)Wn(f3x)} 
a2  —  P2 

xpZn(ax)Wn-i(Px)  —  xaZn_i(ax)Wn(Px) 
a2-  p2 


completing  the  proof. 
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Appendix  F.  Fields  Within  a  Coaxial  Transmission  Line 

The  fields  within  a  coaxial  transmission  line  are  well-documented,  [10].  There  specific 
forms  are  developed  here  as  a  reference. 


F.l  TM  Fields  in  a  Coaxial  Line 


TMZ  fields  in  region  I  and  III  of  figure  3.1  are  generated  from  a  z  directed 
potential  function  that  is  a  solution  to  the  wave  equation  [3].  The  potential  function 
is  (B.15)  and  the  generation  functions  are 


EP  =  ~j 


1  d2T 

u/jle  dpdz 


HP 


1  1<9T 

p  p  dtp 


Ea,  — 


-J 


i  i  a2T 


ujpe  p  d(j)dz 


H4) 


1  dT 

p  dp 


Ez  =  -j- 


ujpe 


Hz  =  0 


which,  considering  only  a  forward  propagating  wave,  produce 


EP 

E^ 

Ez 

HP 

H(p 

Hz 


j—y'm  (kPnp)  [Cm)n  cos  +  Dm,n  sin  e  lz 

Vfl'l  (  1  \ 

j -  -  Ym  ( kpnp )  [Cm>n  sin  (m(f>)  -  Dm,n  cos  (m^)]  e 

\pj 

k2 

( kPnp )  [Cmtn  cos  (m<j>)  +  Dm>n  sin  (m<f>)\  e^z 

UJ  /i€ 

TIT  (  1  \ 

—  -  Vm  ( kPnp )  [Cm,n  sin  (m<j>)  -  Dm  n  cos  (m4>)\  e~lz 

P  \PJ 


-1  z 


k 


— —F'm  (kpnP)  [ Cm,n  cos  ( mq !>)  +  Dm,n  sin  {m4>)\  e 


-7* 


p 


103 


where 


Ym  (kpnp)  lA-m,nJm  (kp„P)  "h  BmriYm  {kpnP)\ 

{kPnp)  —  [ ArnnJm  ( kPnp )  +  BmnYm  {kPnp)\ 

and  m,  n  identify  the  TMmn  mode.  Boundary  conditions  require  the  tangential  electric 
fields  to  be  zero  at  the  inner  and  outer  conductor  of  a  coaxial  line.  Assuming  the 
conductor  is  a  PEC,  knowing  that  the  condition  must  be  met  independently  of  0, 
and  that  the  surface  normal  reference  for  the  tangential  fields  is  p,  the  appropriate 
boundary  conditions  are 


E</,(p  =  a,  0,  z)  =  E^p  =  b,  <j>,  z )  =  Ez(p  =  a,  0,  z)  =  Ez(p  =  b,  <j>,z)  =  0 


where  a  and  b  are  the  radius  of  the  outer  surface  of  the  inner  conductor  and  the 
inner  surface  of  the  outer  conductor  respectively.  Applying  the  boundary  condtion  at 
p  =  a,  solving  for  the  constant 


A  —  —  R 

1  Lm..r.  JSr 


Ym  (kpna) 

,n  Jm  (kPna) 


and  substituting  modifies  the  field  equations  to 


EP 
E # 
Ez 

HP 

H« 

Hz 


r)kpBmn  . 


uj  jj  p  J ni  ( kPn  a  } 


Y'm  (kpnP)  [Cm,n  cos  (mcj))  +  Dm>n  sin  (mcj))]  e 


~lz 


J 


m^Br 


(kPna)  \p 


~  Ym  ( kpnp )  [Cm,n  sin  (mcj))  -  Dm>n  cos  {m4>)\  e 


~lz 


k 2  B 

~j _ Pr  7,'n  ,/rn  ( kPnp )  [Cm,n  cos  (m<f>)  +  Dmjn  sin  (m0)]  e_7z 


L ( kpnd 


mB, 


m,n 


pJm  (kPna)  \p 


~  Ym  ( kPnp )  sin  {m<f>)  -  £>m,n  cos  e 


—Jz 


kn„  B 


Pn  771,77.  t  r/ 


pJm  (kPna) 


Y'm  ( kPnp )  [Cm>n  cos  (mcj))  +  Dmjn  sin  e  7z 


0 
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where 


(kpnP)  [Jm  (kpnp)  (kpn(Y)  ^ rn  (kpnCl)  Ym  (kpnp)] 

Vm  (kPnp)  —  [Jm  ( kPnp )  Ym  (kPna)  —  Jm  ( kPna )  Ym  ( kPnp)\ 

Applying  boundary  conditions  at  p  =  b  produces  the  characteristic  equation  necessary 
to  solve  for  kPn  as 

[Jm  ( kPnb )  Ym  ( kPna )  -  Jm  ( kPna )  Ym  ( kPnb )]  =  0  (F.l) 

where  the  subscript  n  denotes  the  nth  zero  of  (F.l). 

F.2  TEM  Fields  in  a  Coaxial  Line 

TEM2  fields  are  the  dominant  mode  in  a  coaxial  line  and  have  a  cutoff  frequency 
equal  to  0  Hz.  The  development  of  appendix  C  can  be  continued  to  define  the  electric 
and  magnetic  field  equations,  but  a  more  interesting  method  is  to  apply  the  result 
E  —  'Vt'ip  and  kp  =  0  to  the  TM2  fields  derived  in  the  previous  section.  Physically,  a 
potential  gradient  will  not  exist  in  the  (f)  direction  for  a  TEM  field  because  the  con¬ 
ductors  are  coaxial  and  therefore  E  =  implies  a  p  directed  electric  field  potential 
gradient.  Additionally,  the  electric  potential  gradient,  Vt^,  must  be  independent  of 
(f)  making  the  resulting  field  axially  symmetric.  The  TM2  fields  are  axially  symmetric 
when  m  =  0. 

Taking  the  limit  as  kPn  — >  0  for  axially  symmetric  ( m  =  0)  TM2  fields,  using 
the  small  argument  approximations  shown  in  table  F.l,  and  considering  only  forward 
propagating  waves 


E,], 

=  0 

m  =  0  for  axially  symmetric  fields 

HP 

=  0 

m  =  0  for  axially  symmetric  fields 

Hz 

=  0 

TMZ 
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lim  E„  = 


jupe 
7 

jvpe 


kPnp\  /2  ln  hkPna 


(1)l 


(  -2 


^PnP 


o-lnZ 


2kPnp 

nap 


lkPna  ln  ( 7 kPna 


2 

+  [  — 

7 ipy 


?-lnZ 


2kBo,n  1 
nupe  [p 


o~jkz 


lim  k2  =  k2pn  +  k2  —  J2  pe 


lim  E?  = 


PAP 

jujpe 


o~lnZ 


kPn  Bq trl 

jupe 


^lkPna 

2nap 


ln 


7  kPna 


47 kPnp  ^  / 7 kPnP\ 
2vrp7  V  2  / 


■,-lnZ 


=  0 


lim  Ha,  = 


kpnB0n 

9 


^kPnp' |  ^2  ^  ^7fep„« 


(1)  te) 


?-lnZ 


BL 


f  jpp 

p  [  V  7TA7 


7^0 


In 


7^’pn« 


2 

+  (  — 
Tip  / 


o-'InZ 


2Bja 

Tip 


o-lnZ 


/  kzu>ep\ 

2  B*n 

'1' 

\kzLvepJ 

Tip 

.9. 

\  u>ep 

'1' 

\  TTLUep  , 

I  puj^fep 

.9. 

o-lnZ 


D-jkz 


lim  k 2  =  k2n  +  k2  =  uj2pe 


/2  kB*n\ 

1 

'1' 

\  Tiioep  ) 

Ztem 

.9. 

?—jkz 


where  ZTEM  —  kJL 
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Table  F.l:  Bessel  and  logarithmic  function  small  argument  limits 


lim  J0  (x) 

x—>0 

~  1 

X 

lim  Y'0  (x)  ~ 

x — >0 

!•»(?) 

-2 

lim  J\(x) 

_ 

lim  Y\  lx)  — 

2 

i->0 

7 TX 

lim  xa  In  x 

x—>0 

=  0 

(a  constant,  3?(q 

V 

o 

In  summary,  axially  symmetric  TMZ  fields  reduce  to  a  TEM  field  defined  by 


where 


EP 

H(t> 


T 


IP. 

T 


D-jkz 


Ztem  IP 


D-jkz 


2fcB+ 

7T  u/ie 


(F.2) 

(F.3) 


(F-4) 
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Appendix  G.  Transverse  Mode  dependent  Vector  Coupling  Integrals 

This  appendix  documents  the  development  of  coupling  integral  solutions.  It  does  not 
provide  the  greatest  simplification,  but  does  give  a  general  solution  and  a  simplified 
solution  for  both  inner  and  outer  air  gaps.  Coaxial  line  regions  I,  II,  and  III  are  as 
indicated  in  figure  G.l. 

G.l  Symbols  Defined 

G.1.1  General. 

Vm  ( ap )  =  Jm  (ap)  Y0  (aa)  -  J0  (aa)  Ym  ( ap ) 

Wm  (ap)  =  Jm  (ap)  Y0  (ab)  -  J0  (ab)  Ym  (ap) 

where  Jm  is  a  Bessel  function  of  the  first  kind  of  order  m  and  Ym  is  a  Bessel  function 
of  the  second  kind  of  order  m  (sometimes  referred  to  as  a  Neumann  function). 

G.l. 2  Regions  I  and  III  ( Freespace ). 

7n  =  Kn  ~  ko  ko  = 

2  _  In  rj-  _  rgnkpn 

jueo  j2k0 

where  kPn  is  the  nth  solution  of 

Vo  =  Jo  ( kPnb )  Yq  (kPna)  -  J0  (kPna)  F0  (kPnb)  =  0 
G.l. 3  Region  II  (Material  with  Air  Gap). 


k2pm  = 

7n  +  k\ 

kp2n 

=  fin2  +  kl 

h  = 

k\ 

Zln  = 

rfn 

in 

jve  i 

jue  2 
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where  kp in  and  kp2n  are  solutions  to 

^  pin  ^  In  Vo  (kpinR)  W  i  ( kp2nR )  —  kp2n  ^2  n  V  (kpinR)  (kpinR)  —  0 


and 


Outer  Air  Gap 

R  =  r2 

ei,/Ti  are  sample  properties 


62,  /i2  are  6q,  Po 


Inner  Air  Gap 

R,  —  V\ 


G,  Pi  are  6o,  po 
62, /i2  are  sample  properties 


_  '^kpinr)n 

j2h 


c  _  n kp2nin  62  VOi  (kp2nR) 

j2k2  6i  Vi  (kpinR) 


jj  T^k pin^n  6i  I  1  (kpinR) 

j2k\  e2  Wi  (kp2nR) 


_  '^kp2nrjn 

j2k2 


G.2  TEM  testing  operator  coupling  to  Region  II  TM  fields 


Bln 


Cm 


R  b 

/  e i  •  en  pdp  +  j  ei  ■  enpdp 

a  R 

R  b 

/  hi  •  hn  pdp  +  /  hi  ■  h„ pdf) 


R 


R 


ZnZm 


ex  •  en  pdp +  /  ei  • enpdp 


ZnZ‘ 


2n 


R 


(G.l) 
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First  Integral: 


R 


j  ei  •  en  pdp 


a 


{V\{kP\np)}  pdp 


R 

C  J  Vi(kpinp)  dp 

a 

^V0(kplnR) 

Kpln 


Second  Integral: 

b 

I  ei  •  en  pdp 

R 


b 

D  J  j“j  {wi(kp2np)}  pdp 

R 


b 

D  J  W1(kp2nP )  dp 

R 

D  {  -r^W0(kp2nP) 

l  Pp2n 

- — Wo(kp2nR) 

Pp2n 


b 

R 


Therefore 


B\n  —  j — WQ(kp2nR)  —  — — Vo(kpinR ) 

& p2n  rvpln 

C In  =  7  7  u - Wo(kp2nR)  -  77, - V0(kpinR) 

^Jn^J2n^jp2n  ^n^ln^pln 

in  general  form.  Using  the  characteristic  equation  (G.l),  the  general  forms  are  sim¬ 
plified  to  for  each  gap  scenario  to: 
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Outer  Gap 


G.3 


Inner  Gap 


Bln  ~ 


Cm  = 


THn_ 

.7  2  A;  i 


vpln 

^p2n 


-  1 


Vo(kpinR) 


j2k\Zn 


1  (k 


vpln 


^2n  \k'p2n 


1 

•^lra 


Vo(^pln-R) 


-Bi„.  — 


Cm  = 


Kin 

j2k2 


vpln 
^ pin 


Wo(kp2nR) 


7T7n 

1 

1 

(  kpm.\2 

j2k2Zn 

N 

to 

3  | 

Zm 

1 

Q. 

Wq  (kp2nR) 


TM  testing  operator  coupling  to  Region  II  TM  fields 


R  b 

B-mn  j  ® n  pdp  ~\~  j  C  m  •  6  n  f)  dp 

a  R 

R 

Cm.n.  ~  ^  /  e  i  ■  e pdp  “l- 


. .  .m  >  1 


ZnZ 


n^ln 


ZnZ‘ 


^1  *  &nP  dp 


. . .  m  >  1 


71^271 


R 


First  Integral: 

R  R 

/  ei  •  en  pdp  =  HC  V\{kpnp)Vx{kp\np)  pdp 


=  H  C 


=  H  C 


pkpm ^  l(,kpnp')Vo('kpmp')  pkpnVo(kpnp)  \  i(kpinp) 


U2  _  1,2 
^ pn  ™ pin 


R 


Rkpmk  i  ( kpnR)  Vi){kpin-B)  RkpnVQ(kpnR)  Vi(kpinR) 


1,2  _  1,2 
™ pn  ^ pin 
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Second  Integral: 


R 


R 


er  •  en  pdp  =  HD  V,  (kpnp)W,  (kp2np)  pdp 


=  H  D 


pfcp2nVi(kpnP)^Vo(kp2nP)  pkpnV(){kpnp^W  I  (  ^‘p2nP) 


1.2  _  l.‘2 

^ pn  n,p2n 


R 


=  -H  D 


Rkp2nVi(kpnR)Wo(kp2nR)  —  RkpnVo(kpnR)Wi(kp2nR ) 


U2  —  h2 

r° pn  r<Jp2n 


Therefore 


Brnn  =  H  R{C 


,  kpinVi(kpnR)  Vq  (/Cpin.R)  kpnVfl(kpnR^)  V\(kp \nH) 


U2  _  U2 
^ pn  ^ pin 

kpnVo(kpnR^  W^l(^p2n-R)  ^  ^  - 

u  U2  _  U2  f  ■■■m>  i 

"'pn  " Jp2n  ) 


Cmn 


H  RlC 


kpinVi  (^kpnR^Vo^kpinU)  kpnVo(kpnlV)  V\(kp\nR) 


~2 
v  pn 


kpln ) 


D 


ZnZln  {k7p 

kp2viVl  {kpn-R^)  {kp2nR^)  k  pnV()(k  pnR)  ^V\{kp2n  R')  1 


^n^2n  {k { 


2 

pn 


kp  2n) 


7 


. .  m  >  1 


in  general  form,  which  simplifies  for  each  gap  scenario  to 
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Outer  Gap  (m  >  1) 


B  ,rm. 


'K'Jnk  pn  Tt^fiikpln 


j2/ci 

R 

7,2 

Kpn 

kpv 

pn 


>2  r  2  )  kpinVi{kpnR) V0(kplnR)  -\- 

—  Kp2n  J 


79  r  2  J  kpnV0(kpnR) Vi(kpinR) 

p2n)  K pn  ~  * p2n  > 


^  -  kl 


^y  _  kpn  7t^fnkp2n  p 

mn=  32k0ZnJ2k2Zln 

1  e2  \ 

7,2  J“~p  ,  (1,2  _  1,2  \  )  kp2nVl{kpnR)Vo(kp2nR )  + 

/t™  "'pin  eHftpn  h‘p2n)  J 

7,2  _  1,2  7,2  _  7,2  )  kpnVQ(kpnR)Vi(kpinR) 

ri'nr,  A/  ™  A,  „0„  / 


pn 


Inner  Gap  (m  >  1) 


vpn  rvp2n  ^  pn  rvp2n 


Bm.n. 


/Kr1nfcpn  ^"Tn^p2n 

j2/c0  j2k2 


R  x 


7,2  _  7,2  r,2  _  7,2  )  kp2vYl  (kpnR)  Wo  (kp2nR)  + 

Kpn  hpln  hpn  p2n  J 


£2 


kpn-kp2n  ^i(k2pn  -  k2pln) 


kpn  Vo(kpnR)Wi(kp2nR) 


_  T^^nkpu  'T^r)nk p2n  p 

mn  =  j2k0ZnJ2k2Z2n  X 

£i 


\^(k2n  -  k2ln ) 


/, 2  _  7,2 
"‘pn  n'p2n 


kp2n^  1  (  kpnR)  {kp2nR) 


7,2  _  7,2  U2  _  L.2 

pn  "'p2n  "'pn  "'pin 


kpn  V0{kpnR)W1{kp2nR) 
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2  a 


2b 


a/  -AAA/V- 
-xAAAA- 


-A/VW-  b1 


c/  AANV* 


a2  >■■■ >an 


b  +  b  + 

U2  ,...,Un 


b2  ■  ■■■>bn 


Figure  G.l:  The  coupling  of  transverse  electric  and  magnetic  mode  vectors  is  eval¬ 
uate  at  z  =  0  and  z  =  L. 
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Appendix  H.  Higher  Order  Mode  Excitation  By  Waveguide 

Discontinuity 

Not  all  discontinuities  excite  higher  order  modes.  Considering  a  simple  media,  higher 
order  mode  excitation  is  determined  by  examining  the  coupling  of  one  mode,  that 
of  the  excitation,  to  the  infinite  number  of  possible  field  modes  within  the  system. 
The  general  modal  method  developed  in  the  text  (section  3.1)  will  be  applied  to 
the  scenarios  shown  in  figures  H.l(a)  and  (b)  to  prove  the  excitation  of  higher  order 
modes.  It  is  important  to  recall  that  applying  a  testing  operator  in  modal  analysis  is 
equivalent  to  determining  the  power  coupling  between  two  modes  provided  the  testing 
operators  are  the  fields  in  one  section  of  the  coaxial  line. 

The  transverse  interface  in  the  guided  structure  of  figure  H.l(a)  will  not  excite 
higher  order  modes.  For  this  scenario,  the  mode  dependent  transverse  vectors  on  both 
sides  of  the  boundary  (z  =  0)  represent  a  coaxial  line  with  the  same  dimensions,  but 
filled  with  different  material.  Therefore,  they  are  defined  as 


en  =  pEp(z  <  0)  +  4>E^,(z  <  0) 
h„  =  pHp(z  <  0)  +  fiH^z  <  0) 
en  =  pEp(z  >  0)  +  (, bE^z  >  0) 
hrt  =  pHp(z  >  0)  +  (frH'pO  >  0) 


where  the  expressions  for  Ep,  E^,  HP)  and  H $  are  defined  in  appendix  F.  Applying 
boundary  conditions  at  the  single  boundary  in  figure  H.l(a)  and  assuming  a  single 
mode  excitation,  the  system  of  equations  is 


N 

®1  +  Rn^n 

n=  1 
N 

hi  +  Y,  Rnhn 

n= 1 


N 

Yj  TnAn 

n=l 

N 

YDK 

n=l 
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The  final  step  is  to  apply  a  testing  operator  via  the  integrals 


I  em  ■  {IdS1  (m  —  1,  2, . . . ,  N)  ...  Operator  for  E-field 

cs 

j  hm  •  {Ids'  (m  —  1,  2, . . . ,  N)  ...  Operator  for  H-field 

cs 

It  is  shown  in  appendix  D  that  unique  modes  (i.e.  unique  solutions  to  the  wave 
equation)  that  are  defined  over  the  same  surface  of  integration,  CS,  are  orthogonal. 
Therefore,  the  integrations  in  this  scenario 


rb  2?r 


6jn  "  6 np  dpd(j) 


rin  0 

rb  2tt 


*  &np  dpd(f) 


rin  0 

rb  2tt 


hm  •  h np  dpdcf) 


rin  0 
rb  2-7T 


hm  •  h np  dpdcf) 


rin  0 


0  (m  7^  n) 


0  (m  n) 


0  (m  ^  n) 


0  (m  ^  n) 


Recognizing  that  only  the  n  =  1  field,  ei  and  hi,  carries  power  into  the  system  then 
power  is  only  coupled  into  modes  where  m  —  n  —  1.  Therefore,  only  a  single  mode 
propagates  through  the  system. 

It  is  highly  likely  that  the  coaxial  line  step  at  z  =  0  in  the  guided  structure  of 
figure  H.l(b)  will  excite  higher  order  modes.  For  this  scenario,  the  mode  dependent 
transverse  vectors  on  both  sides  of  the  boundary  (z  =  0)  represent  the  fields  within 
a  coaxial  line  with  differing  dimensions  and  filled  with  different  material.  Therefore, 
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they  are  defined  as 


=  pEp(z  < 

0)  +  < pE^z 

<0) 

Tin 

< 

P<ra 

=  pHp(z  < 

0)  +  (pH^z 

<  0) 

Tin 

< 

P  <ra 

=  pEp(z  > 

0)  +  (frE^z 

>0) 

T  in 

< 

p  <rb 

h  n 

=  pHp(z  > 

0)  +  fiH^z 

>0) 

Tin 

< 

p  <rb 

where  the  expressions 

for  Ep,  Eh 

Hp , and 

are  defined 

in 

appendix  F. 

The  system  of  equations  has  the  same  form  as  the  previous  example,  but  the 
limits  of  integration  are  now  different.  Applying  the  testing  operator  for  z  <  0 


ra  2ir 


'  6 nP  dpd(j) 


rin  0 

rb  2? r 


h m  ■  h nP  dpd(f) 


rin  0 


0  (m  7^  n) 


0  (m  ^  n) 


but  for  z  >  0,  the  fields  m  ^  n  cannot  rely  on  the  orthogonality  proof  of  appendix  D. 
The  surface1  over  which  the  fields  z  <  0  are  defined  is  different  than  the  fields  z  >  0 
and  therefore  it  is  highly  likely  that  the  fields  are  not  orthogonal  or 


rb  2ir  ra  2ir 

em  ■  enp  dpd(j)  +  /  /  em  •  {  0  }p  dpdcfi  ^  0 


rin  0 

n  27 r 


n  o 

r„  2ir 


rin  0 


hm  •  h np  dpdcj)  +  J  J  hm  •  {  0  }p  dpd(j)  7^  0 

rb  0 


(m  7^  n) 


(m  7^  n) 


1While  reference  is  made  to  the  surface  of  integration  for  the  sake  of  comparing  coupling  integrals, 
the  particular  solutions  to  Bessel  integral  are  defined  between  radial  boundaries  (the  circumferential 
boundaries  are  repeating).  Therefore,  it  is  equivalent  to  say  the  particular  solutions  are  defined 
between  a  different  set  of  boundaries  and  therefore  orthogonality  is  highly  unlikely. 
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Figure  H.l:  Transverse  discontinuities  in  a  coaxial  waveguide  (a)  will  not  excite 

higher  order  modes.  The  step-junction  discontinuity  at  z  =  0  in  (b)  will  couple 
incident  (n  =  1)  energy  into  higher  order  (n  >  1)  modes. 


Returning  to  the  system  of  equations  for  the  modal  method  and  using  this  result 
produces  (considering  only  the  electric  field  and  N  =  2) 


ra  27 r  ra  27 r 

e!  ■  e1p  dpd(f>  +  R\  j  J  ei  ei/9  = 

rin  0  Tin  0 

2n  rb  27T 

T\  j  j  e  i  •  §1/9  dpdcj)  +  T2  j  /ei-  e2p  dpdcj) 

rin  0  Tin  0 


when  applying  the  testing  operator  ei  and 


ra  27 r 

R2  J  J  e 2  ■  e2p  dpd(j) 

rin  0 


Tb  27 r  Tb  27 T 

7j  y  I  e2-  exp  dpdcj)  +  T2  j  j  e2  ■  e2p  dpdcj) 

rin  0  0 


when  applying  the  testing  operator  e2  where  only  non-zero  terms  are  retained.  Phys¬ 
ically,  if  the  initial  excitation,  ei,  couples  power  into  (excites)  modes  en  where  n  >  1, 
that  power  is  coupled  into  the  remaining  modes,  en  and  en  for  n  >  1,  in  the  system. 
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